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Preface 



The study of phisics at the level of the fundamental constituents of the universe 
has represented a challenge for the human mind. The clasification of such objects 
by their properties is very helpful for understanding the interactions among 
them. One of these properties, their mass, remains as one of the most intriguing, 
since, up to date, there is no theory that can explain its origin. The manner 
in which we understand the interactions at the fundamental level is by the 
exchange of mediators in the Standar Model (SM) of Particle Physics. All the 
processes are then described in terms of Feynman diagrams, and in order to 
compare with the experimental results, one must calculate the amplitudes for 
a given event, adding the contributions of all the possible ways in which the 
event can take place. It is true that this model satisfactorily reproduces the 
dynamics of the particles through the symmetries of the mathematical objects 
that describe such particles, but it assumes that they are massless. On the other 
hand, in the different high energy laboratories around the world, masses for the 
hmdamentcl constituents of the universe have been measured, in some cases 
very accurately. This reflects the fact of the breaking of the symmetry of the 
SM. In order to conciliate this discrepancy between theory and experiment, the 
quest for a self consistent theoretical model and its experimental corroboration 
that explains the origin of the masses of the fundamental blocks of the universe 
has begun. 

The most famous as well as popular of these models is the so-called Higgs 
mechanism, which explains that the mass of the particles emerges by their in- 
teraction with the Higgs boson, a phenomcnom that involves the Spontaneous 
Breaking of the Chiral Symmetry in the SM. The search for this boson in particle 
accelerators in the US and in Europe has impulsed the technological develop- 
ment, due to the technical difficulties and low budgets to produce and detect it; 
and the theoretical advance, which offers simpler channels towards its discovery. 
At the same time, it has promoted the increase of the number of students and 
scientists dedicated to this branch of physics. In spite of the eficiency on which 
the SM describes the phenomenology of the fundamental particles, there are 
still some problems. 

The Higgs boson has not been discovered experimentally yet. The SM as- 
sumes that this boson is fundamental, but it does not exist a fundamental scalar 
in nature, at least, it has not been observed yet. On the other hand, when we 
try to extend the SM to include some other forces by constructing a Grand 
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Unification Theory, we face the hierarchy problem, and in order to solve it, one 
must perform an unnatural fine tunning of the Higgs parameters. Supersymme- 
try or SUSY takes care of this problem at the expense of doubling the spectrum 
of particles. Fundamental scalars are still there and SUSY keeps on eluding us. 

As an alternative, models have been propossed in which the exsistence of 
a scalar particle is no longer necessary for the symmetry breaking, since it is 
given in the presence of condensates. Some of these models are Technicolor 
and its Extended versions, Top Condensate and Top Color. Predictions from 
these models are obtained by non perturbative calculations, that even in the 
simplest cases are very hard to perform and they are based in severly incorrect 
assumptions. For example, in the Top Condensate Models in the context of 
the Schwingcr-Dyson Equations (SDE)-which we shortly will talk about-, the 
QDC corrections are only calculated in Landau gauge. If one repeats the same 
excercise in other gauges, one finds that the top quark mass depends upon the 
gauge, which is physically senseless. Trying to solve this situation in complicated 
theories like QCD is a formidable problem. In this thesis we discuss similar 
topics in a simple model, QED in 2+1 dimensions. 

We decided to study the Dynamical Breaking of Chiral Symmetry in the 
context of Schwinger-Dyson equations. These equations are an infinite set of 
integral relations among the Green's functions of a Quantum Field Theory, and 
they provide the analitic structure that such functions must have. Their imple- 
mentation is extended to the High Energy Physics as well as Nuclear Physics, 
therefore, the solutions for these equations are of interest for a wide sector of 
physicists around the world. 

Solutions to the SDE are obtained by truncating the infinite tower of integral 
equations, and the main task for the solution of these equations should be to 
incorporate all of the gauge identities of the theory in order to calculate gauge 
parameter independent physical observablcs. To pursue this objective, analy- 
tical studies, up to where the complexity of the theory allows, are performed. 
Such studies are either based on the use of ansatze for the Greens function or 
in their spectral representation. There are also succesful attempts to describe 
the dynamics of the fundamental particles by discretising the space-time, i. 
e., assumming a lattice where the dynamics takes place. Then, the different 
processes are calculated with the unavoidable use of computers, which the more 
powerful, the more useful in this context. 

Even though there have been attempts to solve the DSE in QCD and in 
QED in a four dimensional space-time, such studies face a difficulty inhcrct to 
four-dimensional theories: Ultraviolet divergences. Added to this, the initial 
assumptions to truncate the infinite tower of SDE have not allowed to clearly 
identify the sources of the gauge dependence of physical observables, nor the role 
which play the gauge invariance constraints in the restauration of the gauge in- 
dependence of these quantities, like the Ward Identities for QED of the Slavnov- 
Taylor Identities for QCD, or some others, like the Landau-Khalatnikov-Fradkin 
transformations in QED, which describe the precise manner in which the Green's 
function vary under a gauge transformation. 

Therefore, we decided to focus or research in Quantum Electrodynamics on 
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a plane, or QED3, where, besides of enjoying of the benefits which theories 
of phenomena occurring on planar surfaces offer, we can identify exactly the 
role of the initial assuptions made to solve the SDE and the gauge invariance 
constraints, since QED3 lacks of ultraviolet divergences. 

The special features of the electromagnetic dynamics on a plane are stu- 
died in Chapter 1, while the derivation of the SDE as well as the Ward-Green- 
Takahashi Identities (WGTI) and the Landau-Khalatnikov-Fradkin (LKF) trans- 
formations is carried out in Chapter 2. For a better understanding of the phe- 
nomenon! of the Dynamical Generation of Fermion Masses, in Chapter 3 we 
solve the SDE for the Fermion Propagator with the only assuption that the 
fcrmions interact among them in the simplest known way, that is, assuming 
that the vertex of the intaraction is only the bare one. This allows us to identify 
the sources of the gauge dependence for the two relevant physical observables, 
the Euclidean mass and the Chiral Condensate. This scenario also allows us to 
study the role of the WGTI in the restauration of the gauge independence of 
the above mentioned physical observables. The conclusion that the imposition 
of the WGTI is a necesary, but not sufficient condition to guarantee the gauge 
independence of the physical observables is translated as the necessity for the 
incorporation of other gauge invariance constraints. The LKF transformations 
are the next ingredient to consider, but their implementation is not as simple 
as the WGTI. In Chapter 4 we describe the manner to implement these trans- 
formations as a requirement that the solutions for the SDE must fulfill. Being 
then only left with the initial assumption about the vertex, in Chapter 5 we 
advocated our attention to remove it by constructing the vertex for the electro- 
magnetic interaction making use of Perturbation Theory as a guide. The main 
reason for this is that in Perturbation Theory, gauge identities (WGTI and LKF 
transformations) and the gauge independence of physical observables are satis- 
fied order by order. Then if our initial assumption for the fermion-boson vertex 
reduces to its Feynman expansion in the weak coupling regime, we stand our 
best chance to obtain the correct gauge behavior in the nonperturbative regime. 
Our construction is the most ambitious of its kind, because, besides considering 
the gauge invariance constraints of QED3, it deals with massive fermions in 
the interaction. Our vertex provides the first insight of the nonperturbative 
interaction, because when we write an effective Ward identity for the transverse 
part of the vertex (i.e., a relation between the transverse vertex and the fermion 
propagator) , we are left with an interaction explicitly independend of the elec- 
tromagnetic coupling. Consequences of this construction as well as the possible 
paths for the implementation of this vertex in similar studies are discussed in 
Chapter 6, where besides we offer our conclusions. We keep the hope that the 
technique we developed can be applied in future to alternative models to the 
SM which help us to discover the origin of masses for the fundamental particles, 
where more reliable calculations are required with improved vertex ansatze. 
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Chapter 1 

Quantum Electrodynamics 
on a Plane 

1.1 Introduction 

Even though we live in a four-dimensional space-time -three of its dimensions 
being spatial and one temporal-, many theories are also interesting when they 
are formulated in other dimensions. Phenomena like the Dynamical Breaking 
of Chiral Symmetry in Quantum Field Theories are more easily understood 
when such theories restric their dynamics to less than four dimensions. Three- 
dimensional models are very useful to understand phenomena that happen on 
planar surfaces, where precisely one counts with two spatial dimensions and 
one temporal. They are also helpful to study the high temperature behavior of 
four-dimensional theories, since these theories have as infinite-temperature limit 
their three-dimensional counterparts ^ 12 ■ -"- n sucn cases, a three-dimensional 
theory is interesting not only as an academic model, but because they have a 
direct physical interpretation. QED3 is one of such models. This is the subject 
of our interest in this thesis and we start from an introduction to it in this 
chapter. 

In QED3, the phenomenom of Dynamical Mass Generation has been studied 
as an alternative mechanism to that of Higgs in order to explain how funda- 
mental particles acquire such property (an excelent review can be found in 
Although the Higgs mechanism offers a solution to this question, it is partial in 
the sense that it relates particles to their masses, but it says nothing about the 
value of these masses. Besides, if the Higgs boson is detected in experiments 
at LEP and/or Tevatron, it can take a while before determining whether it is 
fundamental or not. 

The theory we study has a coupling e 2 with dimensions of mass, which allows 
the dynamically generated mass to be related to this coupling, and it is no 
longer necessary to introduce by hand a cut-off mass scale for this model, since 
e 2 defines this scale. Even more, this theory is superrenormalizable, and lacks 
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of ultraviolet divergences, one of the most serious obstacles to study physical 
phenomena in four-dimensional theories. The order of divergence of QED3 is 



where / is the number of external fermions in a diagram, b the number of 
external photons and n the perturbative order of the diagram. It can be shown 
that the Green's functions become finite in the ultraviolet region. However, 
it does not mean that the theory is free from divergences. In our case, the 
problem is that there are infrared divergent integrals, at least in the massless 
theory. That is why one is forced to start from certain assumptions that remove 
these divergences to solve the Schwinger-Dyson Equations (SDE). 

Therefore, we have a very interesting model, from which we can learn pretty 
much about the Dynamical Generation of Masses and the analytic structure of 
the Fermion Propagator in general jS], but which is mathematically more easily 
treatable than in four dimensions. Results we obtained can serve as a guide to 
more complicated theories, like Grand Unification Theories. 

This theory, besides, posseses a direct physical relevance. QED3 with a 
dinamically generated mass has been propossed as a model for two-dimensional 
superconductivity, specially for the discovery of superconducting quasiplanar 
oxides at high temperature, like L<i2CuOi± y YBa2Cug [Hj- There are other 
models that describe bidemensional superconductivity at high temperature, like 
anyon models, but they are problematic, because they lead to parity violation, 
which is not observed in nature. 

There are also lattice studies on the structure of the Fermion and Photon 
Propagators which allow us to understand the Dynamical Generation of Masses. 
In this kind of studies, the search for the relevant critical exponents is more 
easily carried out that in analytical studies in the continuum. Also, there are 
studies of QED3 at finite temperature, due to the potentially predictive power 
of this model to explain early universe phenomena, and it is also considered as 
the first step towards the comprehension of the hadronic structure, since in its 
unquenched version, i. e., considering vacuum polarization effects, the theory 
exhibits confinement . This can be seen in a heuristical manner by considering 
the classical potential 1 



Eh 



va = 3-f 




n , 



(1.1) 



V(x) — i dxo 



— OO 



d 3 q 
(2^3 



e Hq-x+q x ) g 2 




(1.2) 



where 



(<l) 



1 



(1.3) 



q 2 [l + U(q)} 



1 Although the arguments of Green's functions are the momenta squared, in order to sim- 
plify the notation we will consider that G = G(p). 
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is related to the transverse part of the photon propagator as : 



(1.4) 



Setting 11(g) = 0, one obtains 



e 

V(x) = — lne 2 a; , 



(1.5) 



which is a logarithmically confinement potential, e 2 is obviously the electro- 
magnetic coupling. 

In this chapter we will show the dynamics of the electromagnetically interac- 
ting objects on a plane. We start by explicitly showing the three-dimensional 
equivalent to the Lorentz Group. Next we study the Dirac equation explicitly, 
using 2x2 and 4x4 representations for the Dirac matrices. Then, we study the 
chiral symmetry, that we will be considering throughout the next chapters, and 
the discrete symmetries of QED3. In the end, we shall discuss why we prefer to 
study QED3 in its representation 4 x 4 in this thesis. 

1.2 Lorentz Group 

Lorentz Group can be represented by 3 x 3 matrices, and there are 9 parameters 
to classify them. However, the relation 



fixes 6 of these parameters, leaving only three, corresponding to two boosts, 
each along the corresponding axis, and one rotation. 
Explicit representation for the boosts are : 



9af3 



(1.6) 




(1.7) 



(7 71) \ 
10 
jv 7 J 



Obviously 7 = \f\ — v 2 . To obtain the infinitesimal generators, we perform the 
parametrization 





(1.9) 
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These infinitesimal generators satisfy : 



(1.10) 



^>=o 






—i 





—i 























—i 











—i 









which explicitly correspond to the matrices 

Kj = | -i | , (1.11) 

K 2 = 

A rotation on the plane is described by the matrix 

/ 1 
R(6) = cos(9 sinfl I . (1.12) 



(1.13) 




This matrix has the infinitesimal generator 

dR 



J = - % de 



8=0 




Now, these three generators are in one-to-one correspondence with the parame- 
ters of the Lorentz Group. They obey the following commutation relations : 

[K U K 2 ] = i J 
[J.tfi] = -%K 2 

[J,K 2 ] = +iKl (1.14) 

As it is well known (see, for instance, jSj), pure Lorentz transformations do not 
for a group. Let us observe that J is a Hcrmitian matrix, but neither K\ nor 
Ki are. However, with the definitions 

A\ = \K X , A 2 = ~iK 2 A 3 = J, (1.15) 

relations (|1.14(l can be re-written as : 

[A i} Aj] = ie tjk A kl , (1-16) 

where 

!+l if (ijk) is an even permutation of (123) 
— 1 if (ijk) is an odd permutation of (123) 
else. 

Therefore, we conclude that the Lorentz Group on a plane is isomorphic to the 
SU (2) group. Owing to the reduced number of degrees of freedom on a plane, 
Dirac equation can also be studied in its matrix representation of dimensions 
less than 4. We discuss this matter in the next section. 
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1.3 Dirac Equation and the Lagrangian 

In order to solve the Dirac Equation, let us recall that Paul Adrien Maurice's 
idea (see, for instance [9]) was precisely to write down a Hamiltonian with linear 
dependence on do, and consequently on di, in such a way that its most general 
form was 

Hip = (a ■ P + j3m)ip , (1.17) 
with the requirement that it must fulfill the relativistic relation 

H 2 i]j = (P 2 + m 2 )ip. (1.18) 

Dirac found that a and f3 are no numbers, but matrices with the following 
characteristics : 

• oti, cti y «3 commute with each other. 

. al = q| = a§ = P 2 = 1. 

Dirac equation can be written then as : 

(iP-m)ip = 0, (1.19) 

with $ = 7^(9^, where 

y»=(/3,/3a). (1.20) 
These matrices satisfy the Clifford's algebra : 

{*f,-f} = 2!r. (1.21) 

Once the representation for the 7 matrices has been chosen, one can proceed to 
solve the Dirac equation. 

1.3.1 2x2 Representation 

In QED3, only three 7 matrices are necessary to write down the Dirac Equation. 
The problem is that these matrices should be trated in a special manner when 
the number of dimensions of the space-time is other than four. In the three- 
dimensional basis for the Clifford's algebra is given by the monimials ^Oj : 

1, 7°, 7\ 7 2 , 7°7\ 7°7 2 , 7V, 7°7V • (1-22) 

If Ti represents any element of this set of monomials, it follows that the matrix 
7 = 7°7 1 7 2 satisfies 

[7,1^=0, (1.23) 

Then, by Schur's lemma, 

7 = cl , (1.24) 
where c is a constant. We also have that 

(7) 2 = -1 , (1-25) 
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therefore, 



c = ±i . 



(1.26) 



We can see that we have two inequivalent representations for the Clifford's 
algebra, depending upon the choice of the two signs for 7. In the case of 2 x 2 
matrices, Pauli's a matrices represent the Dirac matrices. We start with the 
metric 

10 

g" v = ( -1 I (1.27) 
0-1 



and proceed to identify 0] : 

7° = 03, 7 1 = ici> 7 2 = i°2 ■ 

The (anti-)commutation relations among these matrices are then 



{7^,7"} - 



These relations imply 



-f-f 



2<r 

fJ,V _ 



fJ 



= 3 



and we have the traceology 

Tr [7^] 
Tr 

Tr [Yl v l p ] 
Tr [7 Al 7 ly 7 P 7 (T ] 



, 

2 (g^gP" - gWn' 1 



Let us notice that the Lorentz indices run from to 2. In our case 



1=13 = <? 3 , 



7 



(3ot\ 



7 = P a 2 = i&2- 



(1.28) 



(1.29) 
(1.30) 



(1.31) 
(1.32) 
(1.33) 
(1.34) 



(1.35) 
(1.36) 
(1.37) 
(1.38) 



(1.39) 



To solve the equation, let us proceed in the standard way, that is, let us consider 
first the particle at rest. The corresponding equation reads 



or. 



do — m 



dip B {t) 
dt 



1 \1 / Mt) 

1 J \ M*) 



(1.40) 
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So, we find the solutions : 

^(t) = e- ,m VA(0), fe(t) = e !mi fe(0), (1.41) 
such that the wavefunction can be written as : 

W = ( 5%® ) • <" 2 > 

with 

»i,fn\ — ( 
If we look for independent solutions of the form 

0i(O) (,')• M0)=(° l ), (1.44) 

we find that the solutions for a particle at rest are 

1 




m = ( ttW ) ■ 



Mt)=[ I )e" mt , E>0, (1.45) 



M*) = ( i ) e ' mt ' E < °- 

For a moving particle, we look for solutions of the form 

${x) = e - lx - p u{p) . (1.46) 

Let us observe that 

d^(x) = d^e- lx - p ]u( P ) = -ip^e- iX - p u( P ) = -i P ^(x) , (1.47) 
therefore, we can write the Dirac equation as follows : 

(VlV - m)u(p) = . (1.48) 

Since 



I^Pn = 7% - 7 • p 



£ -(p 2 + V) 
p 2 - zp 1 -£7 



in matrix form we have 



E — m -(p 2 + ip 1 ) \ / 
p 2 — ip 1 —E — m J \ ub(p) 



(1.49) 



(1.50) 



which leads us to the following system of equations for the spinor componets : 

(E - m)u A (p) - (p 2 + ip 1 )u B (p) = 0, 
(p 2 -ip 1 )u A (p)-(E + m)u B (p) = 0. (1.51) 
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? ■ 1 
p — ip 



E + m 
1 



p 2 -ip 1 I ^ 



To solve this system, we choose u A (p) — 1, such that 

ub ip) = 
=> i/>i(x) = 

This is the positive energy solution. If now we choose ub{p) = 1, 

Ua(p) = 

=> M x ) = 



(1.52) 



2 i • 1 
P + ip 



E — m 



(1.53) 



p~+ip 
E-rn 
1 



-ip-x 



which is the negative energy counterpart. 
1.3.2 4x4 Representation 

In this case, we choose the Weyl or chiral representation for the 7 matrices [§] 



7 



CT 3 
-<T 3 



7 



iui \ 2 _ ( i<?2 

-iai J ' 7 _ ^ -io 2 

(1-54) 

To solve the Dirac equation, once again we start with the equation for the 
particle at rest, which reads 



ia 3 
-i<7 3 



do — m 



1 
1 







-ia 3 j^-mc J \ tpB(t) 

and that leads su to the following system of equations for the spinors 

dtpA 



0, (1.55) 
0, 



ia 3 



dt 



mip A (t) 
~mipB(t) 



(1.56) 



If we decompose the spinors into their components in the form 



Mt) 



we obtain now the following system of equations : 



^Bi(*) 



dt 



= -mip Al (t) 



(1.57) 



(1- 
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at 
at 

dip (t) 



= mip A2 {t) , 
= mip Bl (t) , 
= -rmpB 2 (*) 



whose solutions are 



i>A 2 {t) 

V>s 2 (i) 



= e 
= e v 



mt ^(0) , 
Va 2 (0) , 

V>B 2 (0) . 



(1.59) 



-imt 



Then, the wavefunction tp(t) can be written as 



( e- imt 1> A M \ 
' Va 2 (0) 

lM<>) 

V e-""Vi3 2 (o) y 



-.imt 



We can, as we did before, choose independent solutions such that the 
at rest can be described as 



i 1 \ 

o 

o 

V o J 
/°\ 



1 



i 
o 

V o y 

/ o \ 




V i J 



(1.60) 
particle 

(1.61) 



-imt 



When the particle is moving, the Dirac equation is expressed as : 





(T 3 — ia x p x — icr 2 p 2 — m 




— C73 + io x p x + i<r 2 p 2 — m 



ua{p) 
ub(p) 



from where we obtain the system of equations 

[cr 3 - ia x p x - ia 2 p 2 - m]uA(p) 
[—(73 + io x p x + ia 2 p 2 — m]uB(p) 



0, 
. 



= 0, 
(1.62) 

(1.63) 



Once more, if we decompose the spinors into their components 

UA ± (?) 



ua{p) = 



UA 2 (?) 
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we arrive to the system of equations : 

(E - m)u Al (p) - (p 2 + ip^UAzip) = 0, (1.65) 

(p 2 - ip 1 )u Al (p) - (E + m)u A2 (p) = , 

-(E + m)u Bl (p) - (p 2 + ip 1 )u B2 (p) = 0, 

-(p 2 - ip^UB^p) - (E - m)UB 2 (p) = 0. 

To solve such system, we take advantage of the consequences of choosing a 
particular normalization for the spinors components, as displayed below : 



u A l (?) = 1 => U A2 (p) = 



E+m 



U A2 (P) = 1^U A1 (P) = P - E ± ^, 



U Bl (p) = l=> U B2 (p) 



p —ip 
E-m 



P +ip 
E+m 



• U B . 2 (p) = 1 => U Bl (p) 

Therefore, we have the independent solutions 

/ , 1 A 



ui(p) 



p —ip 
E+m 



V o , 




1 



p —ip 
E-m 



u 2 {p) 



u 4 (p) 



) 



( 

V 

/ 

V 



p —tp 
E+m 
1 







p 2 +ip 1 
E+m 
1 



/ 



(1.66) 



Since the solutions are decoupled, we can merge the corresponding solutions for 
positive and negative energy as 

/ 1 . \ 



u P (p) = 



un(p) 



1 

P'-ip 1 
E+m 
p 2 +ip 
E+m 
1 



(1.67) 



where we understand 
u P (p) 



V 

/ p 2 +ip 1 \ 

E—m 
1 

1 

p 2 -ip 
\ E-m 



u N (p) 



N 
J N 



(1.68) 



In this way we represent the solutions to the Dirac equation. 
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1.4 Discrete Symmetries 

We are now in position to discuss the discrete symmetries of Dirac equation. 
For that purpose, we will follow the work of where we use only the 2x2 
representation for the Dirac matrices. We should first introduce two Dirac fields 
ip and <j) which satisfy respectively 

(il a d a - 7 d 2 - m)V(aj) = 0, (1.69) 
(-ij a d a + 7<9 2 - m)(j)(x) = , 

where the index a runs from to 1 and 7 = j ^ 1 ^ 2 . In order to study the 
charge conjugation operation C, we must write the Dirac equation when the 
field tp interacts with an external magnetic field : 

{ij a 8 a - ~fd 2 - e(n a A a + i 1 A 2 ) - m}i){x) = . (1.70) 

It follows that, under charge conjugation, 

{i-f a d a - 7 9 2 + e( 7 Q A Q + i 1 A 2 ) - m}^ c (x) = . (1.71) 

On the other hand, we have 

{ij a d a - 7 9 2 + e( n a A a + i 1 A 2 ) - m}C^ T (x) = , (1.72) 

where C in the charge conjugation matrix and T denotes the matrix transpose 
operation. Comparing these equation, we obtain that 

4> c {x) =C4> T {x) , (1.73) 

which is the usual charge conjugation operation. 

Parity P in three dimensions corresponds to the inversion of one axis, say x, 
because the inversion of both the axis can be obtained from a rotation of tt of 
the plane. Dirac equation modified by a Parity transformation then reads : 

(z 7 °9o - il l di - 7 9 2 - m)i/j p (t, -x, y) = . (1.74) 

For the field 0, we rewrite the corresponding Dirac equation without performing 
the Parity operation as follows : 

(z 7 °d + vfdi + id 2 - m)<j){t, x,y) = 0. (1.75) 

Multiplying by 1 = 7 ° 7 ° those terms with 7 1 and 7, our equation then reads : 

( i 7 O 9o + i7 1 7 O 7%+77 O 7 O 5 2 -m)0(i,x,y) = 0, (1.76) 
(i 7 °do-i-f j 1 7 di-j ^°d 2 -m) ( j>(t,x,y) = 0, 

(i7 7 5o- i 7°7 1 5i-7V2-m 7 ) 7 V(^^Z/) = 0, 
(z 7 °<9 - i^di - 7 <9 2 - m)j°(j>(t, x,y) = , 



(1.77) 
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such that, under a comparison, 

ij) p (t,-x,y) = j°<f>(t,x,y) 



(1.78) 



It is convenient at this stage to take a closed look at the mass terms of QED3. 
If we write the four-dimensional spinors as 



4> 



V>2 



(1.79) 



a Parity transformation acts on them in the following way : 

1p! -> (Tl1p2 , i>2 ~* Cllpl ■ (1.80) 

Obviously, the ordinary mass term is invariant under this transformation, 

mtptp = m(^4V2 + VWi) > (1-81) 
since Pauli matrices are unitary and hcrmitian. The other mass term 



(1.82) 



is invariant under chiral transformations (jl.89|l . but not under Parity, H1.78(l . 

For the Time Reversal r, we will need to relate somehow the fields tp and 
Dirac equation under a Time Revarsal operation for the field tj) reads : 



(-ij°d + i^di - id 2 - mW{-t, x) = 
Using the same reasoning as before, the field 4> satisfies : 



(— i777°9o — J777 1 9i 
(«77°7<9o + «77 1 7C>i 



7^2 — m)C<j) T (t, x 
7^2 — m)C(j) T (t, x 
7C*2 — m)C(j) T (t, x 
(ijj°d + i77 1 <9i — 77i9 2 + ™l)lC<\) T it, % 
(— i7°9o — J7 1 9i + 7^2 — m)jC(t> T (t, x 

m)^Ccf> T (t, x 

J2 - 

(~i7°7°a + i7°7 1 9i - 7° 7 a 2 - mj°)j°jC4> T (t,x 
(-«7°9o + rfdi - ld 2 - m) 7 °7C^ T (t, x 

from where we conclude that 

r(~t,x)= 1 a 1 C$ r (t,x) . 



(-ij°d - z 7 1 7 7°ai + 77°7°d 2 
{-ij°d + i7°7 1 7°9 1 - 7°77°<9 2 - m)jC0 T (t,x 



(1.83) 



(1.84) 



(1.85) 



We see that if we require the massive spin- 1/2 fields to be invariant under C, 
P and r respectively, we need to introduce the field <fr which satisfies Dirac 
equation with the second representation for 7. It is worth to mention that if 
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originally the lagrangian does not include such field, as in our case, the usual 
mass term breaks P and r, but retains CPt. 
Finally, if we consider the transformation 

<j) -> <j>' = i<y<f> , (1.86) 

we observe that cf>' satisfies the same equation as <f>, except for the sign for the 
mass term. The boost along the y axis and the rotation around it for the field 
<p correspond to the same transformations for the field ip, but in the opposite 
direction, due to the sign of 7. On the other hand, 4>' behaves as ip under 
Lorentz transformations, although with the opposite sign for the mass term. 
This corresponds to the fact that both the representations of the Clifford's 
algebra are the same (under transformations of 7) as the representation of the 
Lorentz Group. Therefore, we can conclude that if we require invariance under 
C, P and r, there must exist two fields whith the opposite sign for their mass 
terms, and under C, P and r, such field are interchanged. 
The invariant lagrangian under discrete transformations is given by 

C = i^(j°d a + i~/d 2 )ip + mifrip + i4>\j a d a + ijd 2 )4>' - mfifi . (1.87) 



1.5 Chiral Symmetry 

As we mentioned before, we can choose the 2x2 representation for the 7 M 
matrices as the Pauli matrices, and we then use two-dimensional spinors. How- 
ever, there is no other 2x2 matrix which anticommutes with the a matrices, 
therefore, the massive theory does not posses a greater symmetry than the 
massles one This is an obstacle to define chirality. That is why we use 

four-dimensional spinors and also those 7^ matrices from the four-dimcsnional 
space-time. The massless theory in this case is invariat under two chiral-like 
tranformations : 

ifj -> e M7 V , (1.88) 

^ ^ . (1.89) 

and, therefore, the lagrangian is invariant under a global U (2) symmetry with 
the generators 

1,7 3 ,7 5 ,[7 3 , 7 5 ] • (1-90) 

This symmetry, however, is broken with a mass term of the form rmpip. Also, 
a dynamically generated mass will break this symmetry, as in four-dimensions. 
In QED3 there is another posible mass term which is invariant under Chiral 
transformations, but not under Parity (as we saw in the last section). Such 
term has the form 

im^[7 3 , 7 5 ]V> . 
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We are only considering the typical mass term, since it has been shown, by 
analyzing the effective potential, that solutions to the SDE with this mass term 
are energetically preferred |1 II IT^] . Another reason to do so is the conservation 
of parity in QED3. As a result, the Lagrangian that we shall be considering in 
this thesis is the one we are familiar with in QED4, i.e., 

C = ${vfdp - m)V - jF^F^ - ±-(d^) 2 . (1.91) 
4 24 

where the notations carry the usual meaning. 

These are the main features of Quantum Electrodynamics on a plane. To 
study the phenomenon of Dynamical Mass Generation, the one we are concern- 
ing with, we must firstly study the gauge structure of QED3, in particular we 
must know the Schwinger-Dyson equations in this context, and also we must 
explore two of the consequences of gauge covariance of the theory: the Ward- 
Green- Takahashi Identities, which relate Green's functions among them, and the 
Landau-Khalatnikov-Fradkin transformations of these functions, which realize 
their gauge behavior under a variation of gauge. As the Lagrangian does not 
change its form, the derivation of the Schwinger-Dyson equations etc. does not 
really differ from the one in 3 spatial dimensions. We take up these derivations 
in the next chapter. 



Chapter 2 

Schwinger-Dyson Equations 
and Gauge Invariance 
Constraints 



2.1 Introduction 

As in any other theory, the electromagnetic dynamics can be obtained, in arbi- 
trary dimensions, from the lagrangian and its corresponding action 

S = J d d xL(<j>{x),d^{x)) . (2.1) 

The equation of motion for the field <f> is obtained after impossing the staticity 
condition on the action 

SS = . (2.2) 
From this condition, we obtain the Euler-Lagrange equations for this field 

In this chapter we start from the QED action. Functional derivatives of this ac- 
tion lead us to the Schwinger-Dyson Equations (SDE) for the Green's functions. 
We also repeat the derivation of the Ward-Green- Takahashi Identity (WGTI) 
and the Landau-Khalatnikov-Fradkin (LKF) transformations, two of the gauge 
identities of QED necessary in order to ensure that solutions to the SDE in the 
study of the Dynamical Generation of Masses with a Non Perturbative Vertex 
reproduce gauge parameter independent physical observables, as they must be. 
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2.2 Electromagnetic Action 

Electromagnetic quantization can be obtained either from the canonical formu- 
lation or by means of the Feynman's Path Integral. We prefer the later scheme, 
since the integral formulation allows us to neatly obtain the SDE. Richard's 
idea was that in order to know the transition amplitude between two quantum 
states, we must sum over all possible histories in which such transition can take 
place, i. e., 

U(t f ,x f <- U,Xi) = J VxVpe lS , (2.4) 

where T>x and T>p are the measures. 

We saw before that the electromagnetic dynamics is obtained from Dirac 
equation (|1.19(l 

(i^d^ - m) ip = . 
This equation can be obtained from the free lagrangian density 

C free = ^( id ^ ~ • ( 2 - 5 ) 

We must consider the interaction between fermions and the electromagnetic 
field, given by the term cy^A^. To make sure that the lagrangian density being 
invariant under global gauge transformations, we must replace the ordinary 
derivatives by covarian derivatives by using the minimal substitution principle : 

dp ^ Dp = dp- ieA^ . (2.6) 

In this way, the lagrangian density is expressed as : 

tint = <K*7^ ~ • ( 2 -7) 

We still need to consider the interaction of the magnetic field with itself. This 
is obtained from the term 

- \f,^ . (2.8) 

Gathering terms, QED action is 



S$,4,A li ]= / d d 



N 



(2.9) 



where / is a flavor label and iV is the number of different fermion flavors. 
Although we are considering bare quantities, the results we obtain are also valid 
for renormalizad quantities, provided this process is performed properly. As for 
the Green's functions, its bareness will be shown explicitly in order to distinguish 
them from the corresponding complete functions. 

Let us turn our atention to the gererating functional 

Z[fj,r,,J,} = f^,i, t A)e^'*^S^f(^^)^^]) , (2 .10) 
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where ry , ry and J M are the sources for the fermions, antifermions and photons, 
respectively, and where it is defined, as in |3] 

dfi(ip, ip, A) = Il f Vip f VipfU^VAft . (2.11) 

To complete the operational definition of QED, let us note that the action is 
invariant under the abelian local transformations 

f/>(x) -» ip x (x) = e~ ieX{x) ip , 

${x) -» ip x (x) = e ieX(x) ip (2.12) 
A^x)^A x (x) = Ap(x) — dfj,X(x) 

where X(x) is an arbitrary scalar function and, for the time being, we have 
supressed flavor labels. Under such circumstances, the generating functional is 
senseless, since for every single fields configuration {ip(x), ip(x), A^x)}, due to 
the gauge inavriance, there exists an infinite number of related configurations 
{ip x {x), ip (as), A x (x)}, which have the same action 

S$^,A li ] = S$ x ,tl> x > A x ] . (2.13) 

The Grassman integration over ip and ip yield the same result, independently 
of X(x), since the corresponding Jacobian is unity. Then, there is a divergence 
in the functional integration over the field A^. The correct definition of the 
measure must ensure that the integration over tha gauge field is extended only 
to inequivalent configurations under gauge transformations. 

This problem can be solved by introducing the Fadeev-Popov determinant. 
The net effect of this procedure in QED is simply to introduce a gauge fixing 
term in the action. A commonly used choice for this term is 

S$,1> t AJ^SS,1>,Av]=S$,il>,A (l \-±J d d x(d^)\ (2.14) 

where £ is the gauge fixing parameter. 



2.3 Schwinger-Dyson Equations (SDE) 

It is known, from some time ago, that from the field equations for a Quantum 
Field Theory, it can be derived a system of coupled integral equations which 
relates the Green's functions of such theory among them. This infinite tower 
of equations is known as the Schwinger-Dyson Equations [T3|. We will use the 
integral functional formulation to derive the SDE following the works on |3] 
and [Hj. 

We start from eq. I|2.14|l and the gererating functional (|2.1(J|) . Let us take 
into account that the fermionic fields {ip, ip} and their sources {fj, rj} are ele- 
ments if the Grassman algebra, that is, all of these field anticommute among 
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themselves; and that and its source J M are c-numbers. We also take the 
standard notations and conventions, where 

4 = A li -f= S rA lt >y u1 {-f,-f} = 2gT t etc. 

Let us note that for an electron, the physical charge must be e phys — — e, where, 
by definition, e = |e| is the magnitude of the charge of the electron. 



2.3.1 SDE for the Photon Propagator 

Let us consider the generating functional l|2.10|l . The generating functional for 
connected Green's functions ^ [77 , 77 , J^] is given by 



Z[fj,r), J„] = e g ^' J ^ . 



(2.15) 



To obtain the SDE corresponding to the Photon Propagator we simply use the 
fact that the functional integral of a total functional derivative vanishes with 
the apropriate boundary conditions. For example, 







SA^x) 



SS e 



SA^x) 



(2.10) 



SS e 



SA^x) 



i5n iSrj ' iSJ^ 



+ J^x)\Z%n,J tl ] . 



Differentiating the action (|2.14(l . we immediately obtain 



5S ( 



SAJx) 



d P d p g^ 



1 - - );)"<>'■ 



A v + e / -0 / 7 A '-0 / 



(2.17) 



from where it follows that, after we divide by Z ', we can write eq. (I2.17|) as 



dJdftT - ( 1 - 7 I d»d" 



sg 



i5J v (x) 



V f ( — 



SG 



+ 



f 



\5r)f (x) Srjf (x) Sr/f (x) 



sg 



Srjf (x) 



-J*(x) .(2.18) 



This equation represents a compact form of the non perturbative equivalent to 
the Maxwell's equations. This is useful for us to obtain an expression for the 
photon vacuum polarization. Now we can take the Legendre transformation and 
introduce the generating functional for one-particle irreducible (1PI) Green's 
functions, T\tp, ip, A^] : 



G[fj, 7], Ap] = iTty, Ap] + i / d d x [V>V + + 



(2.19) 
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From the Grassman integration, it follows that Z[^>, ip, and consequently 
G[fj,t], Jfj] depend only on even powers of fj and 77, which in turn implies that 
setting fj = rj = after taking the derivative of Q (or Z), we will have nonvan- 
ishing results only for the same number of derivatives with respect to fj and 77. 
Similarily, in the absence of derivatives with respect to the fermionic fields, it 
can be seen that only an even number of derivatives of Z and Q with respect to 
J M survive when we set J = 0. From the last expresion, we have 



J M (x) 



SG 



i5J^(x) 
ST 



SA»(x) 



7/ (x) = 



SG 



iSrf (x) 
ST 



6i/>f(x) 



ip * (x) = - 
fj* (x) — 



SG 



iSrjf (x) 
ST 



(2.20) 



x 



From here we obtain expressions for ip, ip and in terms of 77, 77 and J M and 
viceversa; for example, 



¥ a ( x ) = ¥ a [v,V, J^] = i 



SG[y,y, J M 
Srjl 



with the spinorial indices explicitly shown. It is now easy to see that, setting 
J = after we differentiate T, we will have nonvanishing results only when we 
have the same number of derivatives of ip and ip, in analogy with the case of G- 
Making use of the expresions l|2.21[l . let us take a look at the following term : 



d d z 



S 2 G 



5 2 T 



rj—rj—tp—ip—O 



^ S^(z) Srf (y) 
Z 6 V l(x) 



5-q f a {x) 



7}—rj—ip—'ilj—0 

5 a /3Sf g S d (x - y) 



(2.21) 



i/>=i/>=0 



Therefore, when the fermionic sources (77, 77) are null, we can write eq. I|2.18|l 

as : 



ST 



SA»(x) 



Wh,,, - ( 1 - J 



d^dv 



A»{x) 



~i e^Tr[j^Sp(x, x, [A^])] , 
f 



where we have identified the term 
S f F {x,y,[AS) 



sg 



SG 



Srjf(y)Sr)f(x) 5r)(x)8rf (y) 



(2.22) 



(2.23) 



as the Fermion Propagator of flavor / in an extermal magnetic field A^. One 
of the consequences of eq. (|2.21H is that the inverse of this Green's function is 
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given by 



8 2 T 



8if>f(x)6i/>f(y) 



(2.24) 



^>=i/>=0 



Obviously, the complete Green's function for the fermion Sp(x,y) is obtained 
after setting = in eq. I|2.24|l . 

To obtain the corresponding SDE for the photon polarization tensor, we only 
need to act with S/8A„(y) on eq. 12.22fl and set J M (x) = 0. Let us see : 



<5 2 r 



8A»(x)8A v {y) 



dpd p g^ v - [ 1 - - ) d^dv 



-i e^Tr 
f 



8 d (x - y) 
8 2 T 



7m 



8A v {y) \ 8tpf(x)8^f(x) 



.(2.25) 



The right hand side of this equation can be interpreted in a better way by 
observing that : 



s 2 r 



8A v {y) \ Sipf (x)8ipf (x) 



w=U'=0, 



d d u d d w 



8 2 T 



8ipf (x)8ipf (w) 



8 2 T 



</'=</'=o y 
S 2 T 



8A u {y) 8ipf(u)8ipf(w) \ Sipf (w)8ipf (x) 
an analogous result to 
d 



(2.26) 



ii' = V'=0y 



dx 



[A(x)A- 1 (x) ^ I}= 



dA(x) ._ , .. .dA- 1 (x) 
'A 1 (x) + A(x)- 



dx 
dA~Hx) 
dx 



dx 

-A-\x)^A~\x),{2.21) 
dx 



which hold for finite-dimensional matrices. Equation 12.26fl involves the fermion- 
boson Vertex of [1-PI] 



frf( \ - 6 5 T 



(2.28) 



which should not be confused with the generating functional T. Similarily to 
the fermionic case, the second derivative of T with respect to A^ generates the 
inverse of the photon propagator (A -1 )^ (x, y) (left hand side of eq. $2.2b\ ). 
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Therefore, from eqs. (|2.22l) . (|2.25|) and (|2.26|) . we obtain the SDE for the inverse 
of the photon propagator : 



<5 2 r 



SA^(x)8A v (y) 
1 



1 



d 



d^dv 



5 a {x-y)+U^{x,y) , (2.29) 



where we have identified the photon polarization tensor 11^^ : 

n M „(x, y) = i^^(e f ) 2 / d d z 1 d d z 2 Tr[-f t _ l S f F (x,z 1 )Tl(y; z ll z 2 )S f F {z 2l x)} . 
f J 

(2.30) 

Making use of the traslational invariance, we can write the photon propagator 
in momentum space 



A M „(g) = 



+ {<ln<lv/(q 2 + »e)) 



1 



q 2 + ie 



(2.31) 



1 + U(q) * {q 2 + ^e) 2 ' 

where, as usual, we define the scalar polarization U(q) as : 

H^{q) = {-g^q 2 + q f _ l q^)U(q) . 

Let us note that n(q) is independent of the gauge parametre £ in QED, as a 
result of current conservation. At the lowest order in Perturbation Theory, we 
have that Ii(q) = 0, as well as 

Tl(y;z 1 ,z 2 ) = 1 „5 d {y~z 1 )5 d {y~z 2 ) y (i0-m f )S f F (x, y) = S d (x-y) . (2.32) 

Once we factorized , there is no explicit flavor dependence for the proper ver- 
tex. We have seen that from the second derivative of the generating functional 
T[ip, ip, Afj] we obtain the photon and fermion propagators, and from the third 
one, the proper vertex of the fermion-boson interaction. Generally higher or- 
der derivatives of T[tjj, tp, A M ] yield the corresponding proper Green's functions, 
where the number and type of derivatives yield the number and type of legs in 
the proper Green's functions. The SDE for the Photon Propagator is shown in 
Diagram (1) . 



WW§AAAA/ 



wvwvwvv 



Diagram (1 ) : SDE for the Photon Propagator. 

1 Diagrams were generated with AXODRAW I li 
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The representation of the SDE in momentum space is immediately obtained by 
taking the Fourier transform of the expression in coordinate space, or, more 
easily, using the usual Feynman rules for the diagrams based in the lowest order 
perturbative contribution to the non perturbative quantities. For example, for 
the photon polarization tensor we obtain 

*Wff) = (-DE(e') 2 / ^ 

Tr[(i 7 „)(iS f F (k))(irl(k, k + q))(iS f F (k + <?))] , (2.33) 
where the factor (—1) arises from the fermion loop as usual. 



2.3.2 SDE for the Fermion Propagator 



Following a similar procedure, we can derive the integral equation for the 
Fermion Propagator starting from 







f 6St 



d/i $, ^ A) J^ e {i(SM,i,,AJ + f cf^VW+^l) } (2 34) 

6tp(x) 
6 6 5 



i5n' idrj' i6J 
t/(x) + [ift— + e^7 M 



+ ?/(x) > Z[rj,ri, Jf, 



i6J^{x) ) i5r\(x) 



Zfar,,^]. (2.35) 



The last line of this expresion in the functional non perturbative equivalent to 
the Dirac equation. 

As before, we act with 8/6r]f(y) on this expresion to obtain 



5 d {x-y)Z%r 1 ,J tL ]\ 



— I ijS — m? + e^7 M 



i5J^{x) 



S f F (x,y; [Afj]) = , (2.36) 



rf—rj—O 



with obvious notation. Now, using eqs. (|2.15l) and (|2.21(1 . we can rewrite 



S d (x - y) 



(id- m f + e f A(x; [J]) + eV" 
V *■ 



SJv(x) 



S F (x,y;[A fl }) = 0, (2.37) 



which defines the non perturbative connected two-points Green's function. 

The electromagnetic potential vanishes in the absence of an external source, 
that is, A^x] [J = 0]) = 0, in such a way that it is only written to exhibit the 
content of the remaining functional derivation for eq. I|2.37|l . which can be done 
exploiting the identity (|2.26(l : 



i5J»(x) 



S f F (x,y; [A J) 
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= d 



, d JA v {z) 8 / s*r 



i5J*{x)8A v {z) \ 6i/jf (x)5ipf (y) 



2' J d d z d d u d d w^±S f F (x, u)T v {u, w; z)S f F (w, y) 
e f d d z d d ud d w iA^ix, z)S f F (x,u)T l/ (u,w; z)S f F (w,y) , (2.38) 



- -J 



where in the last line we take J = 0. From here, in the absence of external 
sources, eq. (J2.37JI is equivalent to the expression : 

6 d {x-y) - (i?)-ml)S F {x,y) 

= -i(e f ) 2 [ d d z d d ud d wA^(x,z)j^S F {x,u)T l/ (u,w,z)S F (w,y) . (2.39) 



It is usual to write the SDE corresponding to the inverse Fermion Propaga- 
tor. Therefore, multiplying by S F ~ 1 {y,y')^ integrating with respect to y and 
relabeling y' = y : 

S f F -\x,y) - (i@-mf)5 d (x-y) 

= -i{eJ) 2 [ <f l zAA'"'(i,z)7 (I 4(3:,M)r,(ii I y;z). (2.40) 



The photon porpagator couples eqs. (|2.4U|) and l|2.30|l . In this way, one observes 
that the equations for the two-points functions couple to each other, and both 
depend on the three-points Green's function This is the first indication of 
a general rule which says that the SDE for an n-points function is coupled to 
others of the same order or lower orders, and to functions of order (n + 1). 
Diagram (2) shows the SDE corresponding to the Fermion Propagator 




Diagram (2) : SDE for the Fermion Propagator. 



2.3.3 SDE for the Fermion-Boson Vertex 

The corresponding equation for the three-point vertex can be obtained in a 
similar fashion. For completeness, we present it in momentum space, where it 
is written more concisely : 

zr/( P » = i^ + ^J-0-(is F (p' + i)) 

x (tr^p' + l, p + l))(iS F (p + l))K^(p + 1,^ + 1,1) . (2.41) 

K is the fermion-antifcrmion scattering kernel. The diagramatic representation 
of this equation is shown in Diagram (3) : 
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Diagram (3) : SDE for the Vertex. 



Clearly, couples to the two-point function for the fermion Sf, and to the 
fermion-antifermion scattering amplitude M, a four-point function, which again 
illustrates the general rule. 

To solve the SDE for the Fermion Propagator, the simplest approximation 
for the Vertex that has been used is the so-called rainbow or ladder approxima- 
tion, which consists in approximate M by iterating the lowest order perturba- 
tive contribution to the kernel K, along with the substitution of the fermionic 
propagators by their bare couterparts, S F ^ (p) = \j\i>— m/]. This and other 
approximation will be discussed below, avoiding flavor labels. 



Quenched Approximation 

In massless QED3 in the quenched approximation which corresponds 

to neglect fermion-loop contributions to the vacuum polarization, that is, to take 



in the Photon Propagator, we face infrared divergences with the ordinary Per- 
turbation Theory. A commonly used remedy for this situation is to soften the 
infrared behavior of the Photon Propagator by including fermion-loop contri- 
butions to the vacuum polarization. At the lowest order for a fermion of mass 
m, this contribution to the Polarization Scalar is : 



In a theory with TV massless fermions, the Polarization Scalar is then : 



2.3.4 Solving the SDE 



n(<?) = o 



(2.42) 




(2.43) 



n(fc) = - , 



(2.44) 



where 



Ne 2 



a = 



8 
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such that the photon propagator behaves as l/q for q 2 — ► 0, that is, the infrared 
divergence has been softened without altering the ultraviolet properties of the 
propagator. 

The quenched approximation in the SDE for the Fermion Propagator corres- 
ponds to Diagram (4) : 



Diagram (4 ) : SDE for the Fermion Propagator (quenched approximation ) 

and at the one-loop level in the ordinary Perturbation Theory there is no dis- 
tiction between the quenched and the unquenched approximations. 

Rainbow or Ladder Approximation 

In the study of the Dynamical Generation of Masses, one commonly used ap- 
proximation is to set |3] 



which is known as the rainbow or ladder approximation. If this approximation 
is added to the quenched one, the SDE for the Fermion Propagator decouples 
for the corresponding equations for the Photon Propagator and for the Vertex. 
The following Diagram : 



Diagram (5) : SDE for the Fermion Propagator (quenched approximation and 

bare vertex) 

describes this conjunction of approximations. 
Beyond the Rainbow Approximation 

One of the problems with the rainbow approximation is the violation of gauge co- 
variance, particularly of the Ward-Green- Takahashi Identity. The correct form 
for the fermion-boson vertex is crucial to restore the gauge covariance of the 
SDE and should be such that the above mentioned identity is fulfilled, among 
other requirements. Ball and Chiu have studied the structure of this vertex 
and have propossed their now famous ansatz, which we will discus afterwards. 
The restoration of the gauge covariance for the physical observables is one of 
the main motivations for the construction of the fermion-boson vertex, which 
we will carry out in Chapter5. 




r"(fc, P ) = 7"> 



(2.45) 
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1/N Expanssion 

Ordinary Perturbation Theory in terms of the coupling seems to break down due 
to the infrared divergences of the Green's functions. These divergences can be 
avoided making use of another expanssion parameter. We can take N fermion 
flavors and expand in 1/N 0. It has been shown that massless QED3 is finite 
order by order in this approximation, and besides, the infrared behavior of the 
Photon Propagator softenes, as we pointed out before. Therefore, we have a 
theory with N massless fermions and we take the large- N limit, keeping Ne 2 
fixed in the weak coupling regime. If we write 




(2.46) 



the new perturbative expanssion for the Photon Propagator is shown in Dia- 
gram (6) : 



Diagram (6) : SDE for the Fermion Propagator (1/N expanssion). 
and after its evaluation, it yield a Polarization Scalar for the Photon 



which leads to a 1/q behavior for the Photon Propagator in the infrared domain. 
The SDE for the Fermion Propagator in this scheme is given by 



There exists a controversy on whether in the study of this equation, see for 
example Pennington et. at [20j and Atkinson et. al. 21 , DCSB takes place 
for arbitrary number of flavours or there exists a critical number of such flavours 
separating the chirally symmetric and asymmetric phases of unquenched QED. 
However, we do not take up these matters in this thesis. 

Gauge Technique 

There exist a scheme to solve the SDE which differs substantially in the method 
with the previously mentioned studies : The Gauge Technique [23 ED 123 EH 
1291 1301 HITl 1321 133| . This scheme, based in Minkowski space, assumes that the 
elements of the SDE (Propagators and Vertices) have a spectral representation, 
in term of which the SDE are reformulated and directly solved. For exam- 
ple, it assumes that there exists a spectral function pw, such that the Fermion 




n(<?) = q , 



(2.47) 




(2.48) 
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Propagator can be written as : 

S F ( P ) = dco^l , (2.49) 

and the fcrmion-boson Vertex has a similar form. In fact, an ansatz in this 
scheme is 

/OO 11 
dwp^H-j 7m -j . . ( 2 -50) 
-oo P-U 4-UJ 

Inserting these expressions into the SDE, we obtain a linear equation for the 
spectral density. This is an important feature of the Gauge Technique : It 
reduces the SDE to linear equations. Another advantage is that the Ward- 
Green- Takahashi identity is automatically taken into account. 



2.4 Ward-Green- Takahashi Identity 

One of the consequences of gauge covariance is that Green's functions obey 
certain identites which relate one of these functions to the others. Thes are 
called Ward-Geen- Takahashi identities (WGTI), gl E3 |3J], and they come 
out from the Becci-Rouet-Stora-Tyutin (BRST) symmetry. They play a crucial 
role in the proof for the renormalizability of the theory. One of them, simply 
known as the WGTI, relates the [TPI] Vertex to the propagators, and it has 
been widely implemented in SDE studies based either on the Gauge Technique, 
for intance.[251 H3 |27| EH1 ES 123 EH EH, and on making an ansatz for the 
fermion-boson vertex, (1311371 1351 1551 ET1 1321 EE5| . 

We derive this identity following the textbook starting from the gene- 
rating functional l|2.10[l with the action (|2.14(l . which includes the gauge fixing 
term. Let us recall that without such term (and the source terms), the la- 
grangian corresponding to this action is gauge invariant. This makes Z to be 
infinite and spoils the search for the photon propagator. In order to find a finite 
propagator, we are forced to introduce a gauge fixing term (and a ghost term, 
which in the abelian case, we can absorbed into the normalization). This means 
that the lagrangian related to the action (|2.14|l is no longer gauge invariat. The 
physical consequences of the theory, expressed in terms of Green's functions, 
should not depend upon the gauge, in such a way that Z must be gauge inva- 
riant. This is a nontrivial requirement, and leads us to a differential equation 
for Z ', which we will find below. 

Let us take the transformations (|2.13(l to be infinitesimal, that is, 

if> — » ip — ie\(x)tp 

4> -> $ + ie\{x)i> . (2.51) 

Under these transformations, neither the gauge fixing term, nor the source terms 
are gauge invariant, in such a way that the integrand of Z acquires a factor 

e {i f dx[-±(d> i Ap)d p d p \+J'*d l ,\+ie\(ii4>-i>ri)]} M 52) 
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which, being A infinitesimal, can be rewritten as 
1 



1 + i / dx 



{d^A^dpd" + d^J^X - ie(fjrl> - ^) 



A 



(2.53) 



where we have integrated by parts to remove the derivative operator from A. 
Gauge invariance of Z implies that the operator l|2.53|l . when acting on Z, is 
merely the identity. Since A is an arbitrary function, this implies that 



jd p & , (d»A lt ) + d^J^X - ie(fpl) - i>rj) 



Z = 



Substituting the fields by derivatives with respect to their sources, 

1 S 

we find the following functional differential equation for Z : 



i Sfj 



7 1 6 
i on 



i S 
-8 B p 3^ 



d^J u 



e -q— - n— 
or) or] 







(2.54) 



(2.55) 



(2.56) 



Taking the transformation H2.15J) . the last expression can be written as an equa- 
tion for Q : 

^ (^_^)_o. ,, 57) 

where Q — Q[f)^rj, J M ]. Finally, let us turn this expression into an equation for 
the vertex function T, given by the transformation (|2.19|) . Making use of the 
identities Q2.2ip . eq. I|2.57|l becomes : 



ST 



£ SA^x) 5ip{x) 5ip(x) 



ietp- 



ST 



ST 

ietj; „ - , , = 



(2.58) 



Now, taking the functional derivative with respect to ip(xi) and ip{yi), and 
setting ip = = A = 0, the first term vanishes, and therefore, 



-dft 



s 3 r 



x S^{x 1 )Sip{y 1 )SA^{x) 



ieS(x — x\) 



<5 2 r 



-ieS(x - j/i ) 



S 2 T 



5tp(xi)Sip(yi) 



^=•0=0 

. (2.59) 



The left hand side of this equation is the derivative of the [1-PI] fermion-boson 
vertex (|2.28|) . and the next two terms are the inverses of the exact Fermion 
Propagators (|2.23|) . The content of eq. (|2.59|l becomes clear if we expres it 
in momentum space. For such purpose, we define the proper vertex function 
T M (k,p,q) as 

S 3 T 

dxdxxdy^e^-^-^ 



Sip(xi)5ip(y 1 )5At 1 (x) 



ie{2-n) A 5{q-k-p)TJk,p,q) . (2.60) 
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On the other hand, we define the Fermion Propagator in momentum space as 



dx\dy\e 



i(qx 1 -ky 1 ) 



s 2 r 



5ip(xi)6il>(yi) 



= (27r) 4 6(q-p)tSp 1 (p) . (2.61) 



F 

^=•0=0 



Therefore, multiplying eq. H2.59fl by e l ^ qXl kyi px ' and integrating over x, X\ 
and 7/1 , we have 

qL>r^k,p,q) = Sp 1 (k)-Sp 1 (p) . (2.62) 



or, in the limit k — > p, 



0S7 1 



F 



dp^ 



r M (p,p) ■ (2.63) 



WGTI is one of the requirements for the restoration of the gauge covariance of 
the physical observables that have been employed in SDE studies. Its imple- 
mentation is pretty much simple, and it has been widely used. 

There exist also a WGTI for the Photon Propagator, which is given by the 
expression : 

^n""(g) = 0. (2.64) 

This expression is useful, because it allows us to define the Polarization Scalar in 
the traditional way. The fact that IP" is transverse, leads us to the masslessness 
for the photon. We no longer take into account this identity. 



2.5 Landau-Khalatnikov-Fradkin 
Transformations 

In a gauge field theory, Green's functions transfrom in a specific manner un- 
der a variation of gauge. In Quantum Electrodynamics, and in honor to Lev 
Davidovich and his collegues who firs obatained them, these transformations 
carry the name of Landau-Khalatnikov-Fradkin (LKF) transformations, [441 
1451 146) . These were also derived by Johnson and Zumino through functional 
methods. [471 148| 2 . LKF transformations are non perturbative in nature, and 
therefore, they have the potential to play an important role to address the 
problems of gauge invariance which plague the strong coupling SDE studies. 
In general, the rules governning these transformation are far from simple. The 
fact that they better describe their essence in coordinate space, make them even 
more complex. As a result, these transformations have played a less significant 
and practical role than desired in SDE studies. 

We display its derivation below, following the work of Zumino 08]. We start 
by noticing that in Landau gauge, the Photon Propagator can be written as 

A Ml „(aj;0) = 

2 Fukuda, Kubo and Yokoyama have looked for a possible formalism where renormalization 
constants of the wave function are in fact gauge invariant 1491 



9nv 



dfidv 



d 2 



A c (a:) 



(2.65) 
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where A c (x) is the so-called Feynman function 

5(x) 



AJx) = -■ 



d 2 ' 



(2.66) 



In an arbitrary covariant gauge, the Photon Propagator can be parametrized 
by an arbitrary function A^ in the form : 



A^(x; A d ) = A„ v (x; 0) + d^d v A d (x) . 
Let us recall now that the generating functional is expressed as 
Z[fj,r), J M ] = (0\Te l J 'd^+^+ A ^") | ) . 



(2.67) 



(2.68) 



= . 



Just to obtain the LKF transformations, for the moment we are not assuming 
that 

(2.69) 
(2.70) 



Now, under a gauge transformation, 

Zx%n,J„] = z Q [fje ieX ,r)e~ ieX , J**] e ie I dxJ - 

which can be written in a differential form with the expression 



4! -(""•'' + "5T< I z 



(2.71) 



With these definitions, it is easy to verify that the generating functional satisfies 
the following set of differential equations : 



rid, 



■a 



7 M I ie- 



_5_ 

iSr] 



+ m 



_S_ 

iSfj 

+ m 



Z 

-v^z 
-v}z 
+ x]z 



(2.72) 
(2.73) 
(2.74) 
. (2.75) 



The vector operator a M is introduced firstly for the sake of consistency of the 
notation. It satisfies 



-1 



(2.76) 



and a convenient choice for it defines the different gauges as well. For instance, 

a M = d^-d 2 -ie)~ 1 (2.77) 



defines the Landau gauge, and for a Lorentz frame, characterized by a time-like 
unitary vector n M , 

_ djj, + n^(n ■ d) 
9 2 + (n-d) 2 



(2.78) 
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corresponds to the Coulomb gauge. Turning our attention back to the gene- 
rating functional, if an object J 7 , constructed as a functional derivative of the 
generating functional, is gauge invariant in the sense that it does not change 
with the choice of A, then we have that 

f = 0. (2.79) 
In particular, if we take the Fermion Propagator as 

S F (x,y) = ± - ^ Z , (2.80) 
1/1 br\\y)6r\[x) 

we obtain from eq. I|2.70[l or (|2.71() . after setting rj = fj — 0, 

Z A [0,0,J^]5 FA (a ; ,y) = e[ ie(A(a:) - A(! ' ))+i /^ J '* A ]zo[0, 0,J»]S FO (x,y) , (2.81) 
or, 

i-^-[Z[0,0,J^S F (x )y )} = 

[d^{z) - eS(x -z) + eS{y - z)]Z[0, 0, J>*]Sf(x, y) ■ (2.82) 

To obtain a convenient expression for the generating functional which allows us 
to deduce the LKF transformations, let us consider hrst the Fermion Propagator 
in an external magnetic field _B M given by 

S F [x,y;B,} = S F (B) = -J^ZlL- . (2.83) 
7 M (o M - leBfj,) + m 

We write then the vaccum Polarization Scalar, in obvious notation, as : 

11(B) = e -Tr(WSHB)S^(0)) , (2. 8 4) 

By direct verification we have the identity 



[e lfl§F(B)r 'n{B)} = . (2.85) 



^SB^ iSr] i5r] t 
The generating functional can then be written in the following way : 



z[fj, v ,j»] = e^M^J"n 



s 



e [4(J K +a M a p J p )A c (J'*+a'*a p Jp)-i9 p J p A-| d p J p A d d" ^ 86) 

This expression satisfies eqs. (|2.72() to (12.7411 automatically, and ea. l|2.75[) is 

satisfied when we take = 0. We can verify directly this sentence, except in 
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the case of eq. <|2.72|) . where we should observe firstly that, on acting with the 
operator 



d° d, 



6 



in the last exponential of 



io,j u iSJv , 
we obtain a factor 



J» + a p d p J p 



(2.87) 



Those factors with A and A<j do not contibute. We must now move this factor 
to the left of the terms of Ij2.86|l which contain derivatives with respect to J M . 
The net effect of this operation is to replace 



5 S 

by virtue of l|2.85|l . This verifies eq. (|2.72|l . We also take advantage of 



n 



iSJ u 



d p j p = d p j p n 



iSJ u 



(2. 



(2.89) 



and 



S F 



i5J, 



d p J P {z) = 
[d p J p {z) - eS(x -z) + e8{y - z)]S F 

to find the relation 

Z[0,0,J p ]S F [x,y;J p ] = S F 



iSJ,, 



i5J n 



Z[0,0, J p ] , 



(2.90) 



(2.91) 



which involves the propagator l|2.80|l . This expression is very useful to obtain 
in a direct way the LKF transformation for the Fermion Propagator. We sa- 
tart from considering an infinitesimal transformation of the function A^. Such 
variation induces a transfrormation on the generating functional given by 



(2.92) 



SX 



Now, setting r) = X] = 0, the induced change in the generating functional without 
fermionic sources is 



5Z[0,0,J»] = -- I / d»J li (6A d )d p J p Z[0,0,J fl ] 



or, in finite form 



Z'[0,0, J p ] = e~^U dtlJ ^ SAd)dPJ >'Z[0,0,J p ] 



(2.93) 



(2.94) 
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Therefore, for the Fermion Propagator we have 

s{z[o, o, r\s F ) = l -JJ -^(SA d )^(z[o, o, r\s F ) , (2.95) 

which, along with eq. 1(2.82(1 . implies 

S'p(x V J**) = e ie2 ^ A o( x -y)- SA d(0)] + ieJ[SA d (x-z)-8A d (y-z)]d p J p ( Z )dz 

S F (x,y;J») . (2.96) 

For J = we find 

S' F (X, y- 0) = e ie 2 [SA d (x-y)-SA d (0)] SF ^ y . Q) (2 g7) 

Then, for a finite change of A^, the transformation law for the Fermion Propa- 
gator reads : 

S F (x; = S F (x; o) e ^l A ^)-A d ( x )] ^ g8) 

For the Vertex, we have 

B„(z;x,y\A) = B^z;x,y\0)e^ A ^ A ^-y^ 

+S F (x~y;0)e^ A ^- A ^-y^^-[A d (x-z)-A d (z-y)] , (2.99) 

where is the non-amputated vertex, defined in momentum space in terms of 
the amputated vertex as : 

B^k,p) = S F {k)T v {k,p)S F (p)A IMV {q) . (2.100) 
The transformation rule for the partially amputated Vertex 

A M (fc,p) = S F (k)T ll (k,p)S F {p) (2.101) 
follows from eqs. ((2.67(1 . 1)2.98(1 and ((2.99(1 . and it is simply given by : 

A M (z; x, y\A d ) = A M (z; x, y \ 0)e -^(o)-A d ( x -y)] (2 1Q2) 

In the usual covariant way for gauge fixing, the Photon Propagator takes the 
form : 

which is obtained by taking A^ in eq. ((2.67(1 as : 

A„(:r) = -tfeV-7 (2 - 104) 

where e 2 is the dimensionless electromagnetic coupling, and \i is the 't Hooft 
mass scale in Dimensional Regularization. 
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These LKF transformations are ruled by very complex laws. Being written 
in coordinate space adds to their complexity and they have been less used in 
the context of SDE, as compared to the WGTI, particularly in the study of the 
phenomenon of Dynamical Mass Generation. 

To have a better understanding of the role that either the WGTI and the 
LKF transformations play in the restoration of gauge independence for the phy- 
sical observables, it is necessary to first know the phenomenon we are dealing 
with, and those assumptions that simplify the most the SDE, to our knowledge, 
to make use of the bare vertex. This is the scenario we will develop in the next 
chapter. 



Chapter 3 



Gauge Dependence of 
Physical Observables 



In a gauge theory, at the level of physical observables, gauge symmetry reflects as 
the fact that they be independent of the gauge parameter. Perturbation theory 
respects these requirements and besides the Ward-Green- Takahashi identities 
(WGTI) and the Landau-Khalatnikov-Fradkin (LKF) transformations remain 
valid at every level of approximation. However, this has not been achieved in 
general in the non perturbative study of gauge field theories through Schwinger- 
Dyson equations (SDE) carried out so far although significant progress has 
been made. The gauge technique of Salam, Delbourgo and later collabora- 
tors, EH 03 HH 1213 1110] EH], was developed to incorporate the constraint 
imposed by WGTI. However, as pointed out in [33], gauge technique can be- 
come completely reliable only after incorporating transverse Green functions 
with correct analytic and gauge-covariance properties. Another method widely 
used to explore the non-perturbative structure of the SDE is to make an ansatz 
for t he full fermion-boson vertex and then study the gauge dependence of the 
physical observables related to the phenomenon of Dynamical Chiral Symme- 
try Breaking (DCSB). This method has been quite popular in four dimensional 
Quantum Electrodynamics (QED). For example, the vertex ansatz proposed by 
Curtis and Pennington, [37], has been extensively used to study the gauge de- 
pendence of the fermion propagator and the dynamical generation of fermion 
mass in Quenched QED, e.g., |38l 1391 I4TH 141) . Later on, in the work of Bashir 
and Pennington, |421 I43| . an improved vertex which achieves complete gauge 
independence of the critical coupling above which mass is dynamically gener- 
ated is proposed. These methods use the cut-off regularization to study the 
gauge dependence of the physical observables. As the cut-off method in general 
does not respect gauge symmetry, a criticism of these works has been raised 
recently, |5U] 03 [1)2] . They suggest dimensional regularization scheme to study 
t he chirally asymmetric phase of QED so that the possible gauge dependence 
coming from the inappropriate regulator could be filtered out. 
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Three dimensional Quantum Electrodynamics (QED3) provides us with a 
neat laboratory to study DCSB as it is ultraviolet well-behaved and hence the 
source of gauge non-invariance finds its roots only in the simplifying assumptions 
employed and not in the choice of the regulator. Burden and Roberts, [J], 
studied the gauge dependence of the chiral condensate in quenched QED3 and 
proposed a vertex which appreciably reduces this gauge dependence in the range 
— 1 of the covariant gauge parameter £. Unfortunately, the choice of their 
vertex does not transform correctly under the operation of charge conjugation. 
Moreover, the selected range of values for £ is very narrow, close to the vicinity 
of the Landau gauge. In this chapter, we undertake the calculation of the 
Euclidean mass of the fermion (referred to as mass from now onwards) and the 
condensate for a wide range of values of £ in the bare vertex approximation in 
the followig schemes : 

• Setting F(p) = l. 

• Including the equation for F(p), and 

• Making a partial use of the WGTI. 

We have not achieved a complete gauge independence neither for the mass 
nor for the chiral condensate, which is a sing for the necesity of the construction 
and use of the full vertex. 



3.1 The Fermion Propagator 

In quenched QED3, the SDE for the fermion propagator in the Minkowski space 
can be written as : 

—r»(k,p)S F (kWA%( q ) , (3.1) 



where q = k — p, e is the electromagnetic coupling, T"(k,p) is the full fermion- 



photon vertex, S F (p) and (q) are the bare fermion and photon propagators 
defined as 



s F {p) = i/^ A%( q ) = - 9 -^ + (i-O q -^r, (3-2) 

p q A 

and Sf(p) is the full fermion propagator, which we prefer to write in the fol- 
lowing most general form : 

5 '« - T^m ■ 

F(p) is referred to as the wavefunction renormalization and A4(p) as the mass 
function and £ is the usual covariant gauge parameter. 

Eq. I|3.1[l is a matrix equation. It consists of two independent equations, 
which can be decoupled by taking its trace after multiplying it with 1 and p 1 , 
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respectively. Making use of Eqs. 13.2I3.3|) and replacing the full vertex by its 
bare counterpart, these equations can be written as : 



1 



F(p) 

M(p) 
F(p) 



d 6 k- 



F{k) 



2ttV J k 2 + M 2 {k) q 4 
[ -2(k ■ pf + (2 - 0(k 2 + P 2 )k ■ p - 2(1 - 0*V ] , 



«(2 + Q 
2tt 2 



d A k 



F(k)M(k) 1 
k 2 + M 2 {k) q 2 



(3.4) 



where a — e 2 /(47r) as usual. Carrying out angular integration after the Wick 
rotation to the Euclidean space, the above equations acquire the form : 



1 



F(p) 

Mjp) 
F(p) 



1 



2 

irp z 



dk- 



k 2 F(k) 
l + M 2 (k) 



1 



P 



2kp 



In 



k + p 



k — p 



Tip 



A: 2 +7W 2 (fc) 



k + p 



P 



(3.5) 



(3.6) 



A trivial solution to Eq. (|3.6|l is .M(p) = 0, which corresponds to the usual 
perturbative solution. We are interested in a non-trivial solution by solving 
Eqs. 1)3. 5|) and l|3.6|) simultaneously. Such a solution for M(p) is related to 
the mass m and the chiral condensate (ipip)- Assuming a simple analytic con- 
tinuation from Minkowski to the Euclidean space, neglecting the rotation of 
the integration contour, we define m = A4(m). It is true that this is not the 
physical mass for the fermion, and we do not expect it to be exactly gauge in- 
variant. However, since m ~ M(0), we can consider it as an effective mass. At 
most we can expect this Euclidean mass to be approximately gauge invariant, 
in the sense that it is close to the physical mass 0U]. On the other hand, in 
reference |JJ , Burden and Roberts demonstrated that the standard definition of 
the fermion condensate (Eq. (3.9) in jxj) in terms of an integral over the mass 
function is in excellent numerical agreement with the prediction of the operator 
product expansion [S3] which allows us to write {'t/jip) = Ap 2 M.(p)/(2 + £) (in 



units of e ) in the limit when p 



Such an expansion is valid only for 



values of the gauge parameter in the range [0, 1]. As expected, we find that in 
this limit, M(p) falls as l/p 2 so that the condensate does not depend upon the 
momentum variable p. We shall study the gauge dependence of the mass and 
the condensate in the next section. 



3.2 Effect of the Wavefunction Renormalization 

In studying DCSB, it has been a common practice to make the approximation 
F{p) = 1 so that we only have to solve Eq. (|3.6|l . The justification for this 
approximation stems from the fact that perturbatively F(p) = 1 + 0{a^/ir). If 
a is small and we are sufficiently close to the Landau gauge, one would naturally 
expect that F(p) w 1. Although, it has been quite customary to employ this 
approximation, there exist several works which include both the equations. We 
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Figure 3.1: Mass Function Ai(p) in the F(p) = 1 approxmation. 



study the effects of neglecting the wavefunction renormalization quantitatively. 
Fig. (|3.1() depicts the mass function M(p) for F(p) = 1 in various gauges. As 
expected, the mass function is roughly a constant for low values of p and falls 
as 1/p 2 for large values of p. The integration region chosen is from 10~ 3 to 10 3 
and we select 26 points per decade. The mass probes low momentum region of 
this graph, whereas, the condensate is extracted from its asymptotic behaviour. 
Obviou sly, the mass seems to vary in more or less equally spaced steps with the 
variation of the gauge parameter. In order to obtain a quantitative value of the 
mass, we select neighbouring points p a and pb (p a > Pb)> such that M.{p a ) < p a 
and M(pb) > Pb- We then approximate the mass by the following relation : 



M(p b )-M(p a ) 



Pb - Pa 



(m - p a ) + M(p a ) 



(3.7) 



As for the condensate, the figure does not distinguish between the results for 
various gauges. Therefore, we have to look at the numbers explicitly. Table (1) 
shows the value of the condensate for £ ranging from 1 — 5. Momentum p 



is displayed in units of e and the condensate in units of 10 



The point 



p = 1000 was chosen to calculate the condensate. This number seems sufficiently 
large as the 1/p 2 behaviour seems to set in much earlier (p ss 300), as noted also 
in jZj. In Figs. [|3.2JI and (|3.3|l we display the gauge dependence of the chiral 
condensate and the mass for F(p) = 1 in a wide range of values of the gauge 
parameter. The condensate varies heavily with the change of gauge, roughly 
twice per unit change in the value of £. Gauge dependence of the mass is not 
too different either. 



Gauge Dependence of Physical Observables 



47 



U 10 



CONDENSATE IN VARIOUS GAUGES 

(BARE VERTEX WITH F=l ANDW1TIIOUT WGTI) 



Figure 3.2: Condensate < ipip > in the F(p) = 1 approxmation. 



Repeating the exercise by taking both the equations, namely Eqs. (|3.5[) and 
(|3.6|l . into account, we see similar qualitative behaviour of the mass function. It 
is roughly a constant for low values of p and falls as 1 /p 2 for large values of p, 
Fig. H3.4JI . The large p behaviour is also evident from the entries in Table (2). As 
for the wavefunction renormalization, it also is constant for small values of p. As 
p becomes large it goes to 1, Fig. (|3.5fl . In Table (2) we also give a comparison 
with the work of Burden and Roberts, 0. As mentioned earlier, they restrict 
themselves to the close vicinity of the Landau gauge, where our results are in 
excellent agreement. We investigate the gauge dependence of the condensate as 
well as the mass far beyond the Landau gauge. A graphical description can be 
found in Figs. I|3.6|l and (|3.7|l . The following points are important to note : 



• The wavefunction renormalization plays an extremely important role in 
restoring the gauge invariance of the chiral condensate as well as the mass 
of the fermion. Although the qualitative behaviour of the mass function 
in various regimes of momenta remains largely unchanged, whether or not 
we employ the approximation F(p) = 1, quantitative dependence of the 
physical observables mentioned above on the covariant gauge parameter £ 
reduces a great deal by including the wavefunction renormalization. 

• As we move away from the Landau gauge towards large positive values of £, 
the gauge dependence of the condensate as well as the mass keeps dimini- 
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MASS IN VARIOUS GAUGES 

(BARE VERTEX WITH F=l AND WITHOUT WGTI) 



Figure 3.3: Mass in the F{jp) = 1 approxmation. 



shing, without resorting to any sophisticated ansatze for the fermion-boson 
interaction. 



3.3 Effect of the Ward-Green- Takahashi Iden- 
tity 

The bare photon propagator which appears in Eq. I|3.1|) can be split up in 
longitudinal and transverse parts as follows : 



(3.8) 



where 



A„„(g) = -g^/q +^<W<? 



Employing this decomposition, we can rewrite Eq. (|3.1|l as 

d 3 k 



Sf\p) = S%- x (p)-ie 



(2tt) 3 

d 3 k „„„ , „ x ,, q^qy 



: T l/ (k,p)S F (k)Y 



(3.9) 



It is well known that the use of the WGTI, in the equivalent of the last term 
of Eq. (|3.9(l in QED4, filters out a spurious term which is an artifact of using 
the gauge dependent cut-off regulator. Therefore, one is naturally motivated to 
use this decomposition in dimensions other than four. Now employing the bare 
vertex ansatz T fl (k,p) — 7 M , multiplying Eq. (|3.9|l by 1 and p 1 respectively and 
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MASS FUNCTION IN VARIOUSGAUGES jp " 

(BARE VERTEX WITHOUT WGTI) 5=2 




Figure 3.4: Mass Function incliding the equation for the Wavefunction Rcnor- 
malization. 



Wick-rotating to the Euclidean space, we obtain the following equations : 



1 



F(p) 

M(p) 
F(p) 



1 + 



2ttV 



d A k 



F(k) 



1 



k 2 +M 2 (k) q 4 



2(q-p)(q-k) + 

d 3 k 



a 

2^2 



F(p) 
F(k) 



[p 2 (q-k)+M(k)M(p)(q-p)} 



k 2 +M 2 {k) q 2 



2M ( k ) -^pjf) \M(k)(p-q)-M(p) (k-q)} 



On carrying out angular integration, 



1 



F(p) 



1 + 



vp 2 J 



dk 



p i 
2fc ln 



k 2 — p 2 
+M{k)M{p) 



k 2 F(k)/F(p) 
k 2 +M 2 {k) 

k + p 



M( P ) 
F(p) 



dk 



k — p 
1 

k 2 -p 2 ' 

k 2 F{k) 
k 2 +M 2 {k) 



1 



In 



2kp 

2M(k) 
kp 



k- 






k- 


-P 


}} 


ln 




-P 


k - 


-P 



£ (M(k)-M( P ) M(k)+M( P ) 



F(p) 



k 2 — p 2 



2kp 



In 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



As the terms of the type l/(k 2 — p 2 ) are harder to deal with numerically, we 
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F(p) IN VARIOUS GAUGES £} ' 

(BARE VERTEX WITHOUT WGTi) ^=2 



Figure 3.5: Wavefunction Renormalization. 

use the approximation F(p) = 1 to analyze the effect of the WGTI. Under this 
simplification, we only have to solve 



M{p) 



dk 



k 2 



k 2 +M 2 (k) 



2M{k) 
kp 



hi 



P 



'.M{k)-M{p) M{k)+M{p), 
-? < — s — In 



k 2 — p 2 



2kp 



(3.14) 



The mass function obtained on solving Eq. (|3.14|) is depicted in Fig. I|3.8|l . which, 
along with Table (3), reveals that its qualitative behaviour remains unchanged 
both for small and large values of p. In Figs. 1)3. 9|) and l|3.10|l . we compare the 
gauge dependence of the condensate and the mass with and without the usage 
of the WGTI. We find that the gauge dependence of these quantities seems to 
increase by incorporating the said identity. Similar behaviour was observed in 
QED4 by Gusynin et. al. jH3] in studying the gauge dependence of the critical 
coupling above which mass is generated. They carried out a numerical analysis 
of the criticism raised by Dong et. al. [T^] on the work of Atkinson et. al. 

who did not employ the WGTI as suggested in Eq. I|3.9|l . resulting in the 
appearance of a spurious cut-off dependent term 1 . Gusynin et. al. found that if 
one employs the WGTI, the critical coupling is more steeply gauge dependent. 
We find similar behaviour for the mass and the condensate in QED3 in the 
approximation F(p) = 1 in this section. 



1321 



1 To trace back the origin of this error, consult the footnote on page 7680 of the referenc e 
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CONDENSATE IN VARIOUS GAUGES 

(BARE VERTEX WITHOUT WGTI) 



Figure 3.6: Condensate: Comparison between the cases with and without the 
usage of F(p) = 1. 

3.4 Dimensional Regular izat ion Method 

In this section we compare our numerical results with those obtained by em- 
ploying the dimensional regularization scheme, |50l 1511 152] . For simplicity, we 
restrict ourselves only to the Landau gauge without incorporating the WGTI. 
In this case, the equation for the mass function acquires the following form in 
Euclidean space in arbitrary dimensions : 



d d k M{k) 
(2nf k 2 + M 2 (k) q 2 



M(p)=4xa{d-1) I , d , 2 , )A,,^ ( 3 - 15 ) 



where a is a dimensionful coupling except in four dimensions. We define d = 
4 - 2e and 

a = ad /J. 2£ , (3.16) 

ad being dimensionless. We now use the volume element d d k = k d ~ 1 dkdVld-, 
where d£ld is the d-dimensional solid angle defined as 



1=1 

The angle 9d-i varies from to 2tt, whereas all other angles vary from to 
7r. Choosing 9\ to be the angle between k and p, we can easily carry out the 
remaining angular integrations to arrive at : 

, . , , 2(d-l)a f°° Jl2 k d - 2 M(p) r ,„ sin d - 2 6»i 

M(p) = \ ; ' - / dk 2 — — — f{ / dOx = — ^ .(3.17) 

(4^)^r(^i)io k 2 +M 2 (p)J q 2 
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MASS IN VARIOUS GAUGES 

(BARE VERTEX WITHOUT WGT1) 



Figure 3.7: Mass: Comparison between the cases with and without the usage 
of F(p) = 1. 



Using the standard formula, 55 , 



dx- 



o [1 + 2acosx + a 2 ]' 



integration over 9\ yields : 
(3 - 2e)a 



B(a, l/2) 2 Fi(A, A - a + 1/2, M + 1/2; a 2 ) 

lol < 1 , 



M( P ) = 



(4 7 r) 1 - e r(2 - e) 



dk 



2 {k 2 f^M{k) 



1 



Fx l,e;2-e 



P 



k 2 +M 2 (k) 



6(k 2 ~p') + 



1 / fc 2 

-9 2-^1 I 1, e; 2 — e; 

y p z 



9(k 2 -p 2 



(3.18) 



This equation was studied in detail in [SUj in four dimensions, taking e to 
be a small positive number. The factor (/c 2 )~ e in the numerator regulates the 
otherwise divergent behaviour of the integrand for large momenta. As noted in 
|5()| . the hypergeometric function does not play any role in regularization and 
hence can simply be replaced by F(l, 0; 2, z) = 1. In case of three dimensions, 
the hypergeometric function develops a pole for k 2 = p 2 , as is obvious from the 
following identity : 

1 - e 



2 F 1 (l, £ ;2-e;l) = 



(3.19) 



1 - 2e 

As was pointed out earlier, such terms are hard to deal with numerically. 
Due to increasing computational time and memory, we go only up to e = 0.48, 
starting from e = 0.4. To obtain satisfying results, we need to use increasingly 
more points per decade as we approach closer to e — 0.5. For instance, we use 
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MASS FUNCTION IN VARIOUSGAUGES O ' 

(BARE VERTEX WITH F=l ANDWITH WGTI) 5=2 

5-5 - 



Figure 3.8: Mass Function incliding the effect of the WGTI for F(p) = 1. 

100 points per decade for e = 0.46. The problems of ever increasing computa- 
tional time and memory limited us to use 140 points per decade for the case 
of e = 0.48. Despite this large number, we belive that the corresponding result 
falls short of the desired accuracy. As a result, there is a slight rise at the end 
of the flat region of the mass function, and the final descent begins rather late, 
Fig. (|3.11(l . The problematic pole for e = 0.5 in Eq. (|3.19(l corresponds to the 
relatively well-controlled singularity in the following expression 

/ 1 3 ,\ 11 + 2 , 

aJ H 1 '2 ; 2^) = (3 - 20) 
for z — ► 1. A comparison between the mass function obtained from techniques 
based upon the dimensional regularization scheme and the one computed in 
Section 3.3 is also depicted in Fig l|3.11|l . Taking the numerical limitation for 
e = 0.48 into account, we note that as e approaches the value of 0.5, we get closer 
and closer to the result obtained previously, where we work in 3-dimensions to 
start with. 



3.5 Towards the Full Vertex 

WGTI with the bare vertex is valid only if bare fermion propagators are in- 
volved. When we use full propagators, physical observables calculated from 
the corresponding SDE exhibit gauge dependence. Temptation emerges then to 
move towards the full vertex, and the first requirement we must think of is that 
it should satisfy the WGTI. Hence we face the construction of the Ball-Chiu 
vertex detailed displayed below. The full vertex satisfies WGTI 



(3.21) 
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CONDENSATE IN VARIOUS GAUGES 

(BARE VERTEX WITH ANDWITIIOUT WGTI FOR F=l) 



Figure 3.9: Effect of the WGTI on the gauge dependence of the condensate for 
F(p) = 1. 



This relation allows us to decompose the complete vertex into its longitudinal 
T^(k,p) and transverse Tt^(k,p) parts : 



T»(k,p) = T»(k,p) + r»(k, P ) , 

where the transverse part satisfies 

q^(k,p) = y T^(p,p)=0 



(3.22) 



(3.23) 



and therefore, remains undetermined by the WGTI. Following the work of Ball 
and Chiu, we define the longitudinal part of the vertex in terms of the Fermion 
Propagator. In the limit k — > p, WGTI is written as : 



(3.24) 



Substituting the expresion for the full Fermion Propagator, we observe that 



r"(p, P ) 



d 



dpn 

d i>- M(p) 
W» F{p) 

1" 



d 



jt d M(p) 



f( p ) 1 2p ^d P 2 2p dp 2 f{ p ) 



(3.25) 



After we make the substitutions 



F(p) 



F(k) F(p) 
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(BARE VERTEX WITH ANDWITI [OUT WGTI FOR F=l) 



Figure 3.10: Effect of the WGTI on the gauge dependence of the mass for 
F(p) = 1. 



_9 1_ 

dp 2 F{p) 
d M{p) 
dp 2 F{p) 



\{W+p») 



1 



k 2 — p 2 
1 

k 2 — p 2 



1 



1 



F(k) F{p) 
M{k) M(p) 



F(k) F{p) 



(3.26) 



we can define the longitudanal or Ball-Chiu (BC) vertex as 



BC — 



+ 



F{k) F{p) 



1 {}i + fj{k + pY 
2 



(k 2 -p 2 ) 



(k + P y 

(k 2 — p 2 



M(k) M{p) 



F{k) F{p) 



F{k) F(p) 



(3.27) 



Without loss of generality, the transverse vertex can be expressed as : 



r£(fc, P ) = J2n(k 2 ,p 2 , q 2 )Tt(k,p) , 



(3.28) 



with the apropriate {T M } basis. Functions n, T4, T5 and T7 are proportional to 
the mass m, and therefore in massless studies they do not appear. During the 
last few years, a programme has been started towards the construction of vertex 
ansatze which impose constraints on the transverse part of the vertex. Some of 
the most famous attempts are discussed below. 
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Figure 3.11: Mass Function for various values of e. Notice that e = 0.5 corres- 
ponds to 3 dimensions. The result for e = 0.48 could not be achived with the 
desired accuracy due to the computing time and memory increase. 

3.5.1 Curtis-Penington Vertex 

In four dimensions, the question of why cannot we set all of the t's as zero was 
answered by Curtis and Pennington, arguing to Perturbation Theory and the 
multiplicative renormalizability properties of the Fermion Propagator. 
Their vertex (CP) 



CP BC ' 2 



1 1 



F(k) F{p) 



Y{k 2 -p l )-{k + P Y4 

d(k,p) 



(3.29) 



where 



«M - ^->sf + \M>m + MH P )? (3 30) 

comes out from the following assumptions : 

• Transverse vertex must agree with perturbative results at the one-loop 
level in the relevant kinematic regime k » p. 

• Fermion Propagator obtained from this vertex should be multiplicatively 
renormalizable . 

• This Propagator should transform correctly under its LKF transformation. 

• Vertex should have correct charge conjugation properties. 

• Transverse vertex should vanish in Landau gauge. 

• Transverse Vertex should not depend upon the angle between fermion 
momenta. 
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• The only contributing coefficient is tq. 

• Transverse vertex as well as the longitudinal one are writen without ex- 
plicit dependence on the gauge parameter £ 2 . 

Transverse vertex is written in such a way that, in the massless case, their 



T 6 = 7: 



f k 2 + p 



2 



2 (k 2 -p 2 ) 2 



1 



F(k) F{p) 



(3.31) 



is antisymmetric under the interchange of k and p. 

This vertex, however, does not lead to a completelly gauge independent 
critical coupling. In the massless case, it exhibits a kinematic singularity when 
k 2 — > p 2 . In comparison with Perturbation Theory at one-loop, CP vertex 
agrees with perturbative results only in the kinematic regime above mentioned. 
Multiplicative renormalizability for the Fermion Propagator is achieved only 
in the leading term of the power law for the wavefunction rcnormalization. 
Finally, the LKF transformation of this propagator is valid only up to the leading 
logarithmic term. 



3.5.2 Burden-Roberts Vertex 

Burden and Roberts have parametrized a slight modification to the BC vertex 
in the following way 7 : 



1 BR 



'F{p) 
■a) ft 



aft) 



(k+pY 
k 2 — p 2 



k 2 — p 2 
M{k) 



M{p) 



F{k) F(p) 



1 



1 



F(k) F(p) 



(3.32) 



where 



a=-+S 



(3.33) 



With the assumptions that the vertex should satisfy the WGTI and that it 
should be free from kinematical singularities to all orders in Perturbation The- 
ory, their parametrization allows to optimize 6 such that the most gauge inde- 
pendence of the chiral condensate is achieved. They report S — 0.03 as the best 
value. However, we observe that the BR vertex can be written as : 



1 BR 



1 BC 
1 BC 



1 



1 



F(k) 
St 6 T» , 



F(p) 



[k 2 



{ 7 *V -k 2 ) + (k+p)»i 



(3.34) 



2 This, as we will show, is not possible. 
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with 



T 6 = (fc 2 



V 



1 



1 



(3.35) 



F(k) F(p)_ 

Charge conjugation symmetry for the vertex requires t 6 to be antisymmetric 
under the interchage of k and p, but the BR vertex is symmetric. Besides, for 
S = 0.03 this vertex does not agree with perturbative results, and therefore we 
neglect it as a good choice. 



3.5.3 Dong-Munczek-Roberts Vertex 

Dong et. al. [H] have proposed a vertex (DMR) that satisfies WGTI and 
ensures gauge covariance of the Fermion Propagator in the massless case by 
means of the so-called Transversality Condition : 



(3.36) 



DMR vertex is constructed under the assumption that the coefficients of the 
transverse part are independent of the angle between fermion momenta. In it, 
the following constraints are set on the transverse coefficients : 



for i± 3, 8 
1 (d-l)I 3 



1 h 



d is the number of space-time dimensions and 

h = i((*»+p?)zi-i) 
h = \[k 2 + P 2 ]h, 

where they define 

being the solid angle 
1 



In — I d^d 



{k-pf 



dTld = 



and 



dO-2 sin d " 2 I 



TV 



d6 3 sin d - 3 6» 3 . 



r(d/2) 



2jt 



d0, 



d-1 



(3.37) 



(3.38) 



(3.39) 



(3.40) 



(3.41) 
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Functions p are related to the in the form : 

1 

k 2 — p 2 

In this construction it is shown that the effect of including only t§ as in the CP 
vertex can also be reproduced by considering T3 and t%. This vertex, however, 
exhibit logarithmic kinematic singularities. It has the correct charge conjugation 
properties and leads to a multiplicatively renormalizable Fermion Propagator. 

3.5.4 Burden- Tjiang Vertex 

Using a similar reasoning, Burden and Tjiang |5fi| have deconstructed a one- 
parameter ansatze family for massless QED3. Constraints on the transverse 
coefficients in the BT vertex are the following : 



(3.43) 

\{k 2 +v 2 ) 

\{k 2 +p 2 ), (3.44) 
(3.45) 

Choosing (3 = 1 leads to the DMR vertex. This deconstruction is based upon 
the assumptions that the transverse vertex vanishes in Landau gauge and takes 
tq = 0. Transversality Condition is the strongest argument to ensure the gauge 
covariance of the Fermion Propagator. A possible drawback mentioned by the 
authors is that for any value of (3, for k 2 = p 2 but fc M ^ p^, the BT vertex 
exhibits a logarithmic singularity. 

A criticism to ths vertex has been exposed by Bashir et. al. |F7]. They as- 
sure that T6 cannot be taken as zero and that the transverse coefficients should 
depend upon the angle between fermion momenta. One important observation 
is that the parameter /3, introduced explicitly gauge parameter independent by 
the authors, should have such dependence as pointed out by the corresponding 
perturbative calculation. Finally, in this work they point out that the Transver- 
sality Condition is not valid beyond the one-loop level. 



F(k) F{p) 



(3.42) 



where 



3 i(k, P ) 

1 ' J(k, P ) 

f = 0, 



/(m = (fc2 :/ )2 h/( fc+ ^ 



16kp \ (k — p) 2 

j(M = ( fc2 7^ 2 h/( fc + ^ 2 



IQkp \ (k — p) 2 
They have used the relation (|3.42|) and they define 

f = r 8 + (fc 2 + p 2 )r 2 . 



GO 
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3.5.5 Bashir-Pennington Vertex 

With the assumptions that the transverse vertex vanishes in Landau gauge and 
thet it has no dependence upon the angle between fermion momenta, Bashir 
and Pennington 021 ■ have constructed a vertex that satisfies : 

• WGTI, 

• that the Fermion Propagator is multiplicatively renormalizable, 

• that with the correct choice of the functions that define this vertex, it 
agrees with Perturbation Theory in the weak coupling regime is war- 
rantied, 

• offers a strictly gauge independent critical coupling. 

BP Vertex imposes integral and diferential contraints on the functions that 
define its transverse part. For the massless case, Kizilersii, Bashir and Penning- 
ton |58j removed the initial assumptions, including all the corrections in the 
exponent of the power law for the wavefunction renormalization. This is done 
by introducing new functions and constraints on them. 

Continuing with the programme, we should start an ambitious quest towards 
the construction of a vertex for the fermion-boson interaction in massive QED3, 
taking into account all the nice features of the above mentioned vertices, but 
getting rid of their drawbacks. It is clear that the inclusion of the WGTI is 
common in all the vertex ansatze. This is due to the fact that its implementation 
is straghtforward. Demanding correct gauge behavior of Green's functions is a 
nontrivial requirement. We will show thos in the next chapter, where we study 
the LKF transformation for the Fermion Propagator at the tree level. 
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Appendix: Tables 

Below we display the Tables that contain the numeric results of our calculation : 





V 






P 


^p z M(p) 


0.0 


1000 


2.31109 


2.5 


1000 


11.6992 




642.233 


2.31117 




642.233 


11.7001 




316.228 


2.31119 




316.228 


11.7003 


0.5 


1000 


3.61103 


3.0 


1000 


14.4432 




642.233 


3.61119 




642.233 


14.4444 




316.228 


3.61124 




316.228 


14.4448 


1.0 


1000 


5.19982 


3.5 


1000 


17.4761 




642.233 


5.20009 




642.233 


17.4777 




316.228 


5.20017 




316.228 


17.4782 


1.2 


1000 


5.91622 


4.0 


1000 


20.7977 




642.233 


5.91654 




642.233 


20.7998 




316.228 


5.91664 




316.228 


20.8005 


1.5 


1000 


7.07745 


4.5 


1000 


24.4081 




642.233 


7.07787 




642.233 


24.4108 




316.228 


7.07800 




316.228 


24.4117 


2.0 


1000 


9.24390 


5.0 


1000 


28.3073 




642.233 


9.24452 




642.233 


28.3106 




316.228 


9.24473 




316.228 


28.3117 



TABLE 1. Condensate in Various Gauges for F(p) = 1 
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P 




p 2 M(p) 






BHR 


BR 


0.0 


1000 


2.31109 


2.316 




642.233 


2.31117 






316.228 


2.31119 




0.5 


1000 


1.77309 


1.775 




642.233 


1.77313 






316.228 


1.77306 




1.0 


1000 


1.44791 


1.447 




642.233 


1.44793 






316.228 


1.44780 




1.2 


1000 


1.35288 


1.352 




642.233 


1.35288 






316.228 


1.35274 




1.5 


1000 


1.23591 






642.233 


1.23591 






316.228 


1.23574 




2.0 


1000 


1.09014 






642.233 


1.09012 






316.228 


1.08992 
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P 


BHR 


BR 


2.5 


1000 
642.233 
316.228 


0.98597 
0.98594 
0.98571 




3.0 


1000 
642.233 
316.228 


0.90941 
0.90938 
0.90912 




3.5 


1000 
642.233 
316.228 


0.85201 
0.85197 
0.85169 




4.0 


1000 
642.233 
316.228 


0.80839 
0.80833 
0.80803 




4.5 


1000 
642.233 
316.228 


0.77497 
0.77491 
0.77458 




5.0 


1000 
642.233 
316.228 


0.74933 
0.74927 
0.74891 





TABLE 2. Condensate In Various Gauges Including F(p) 



Gauge Dependence of Physical Observablcs 





V 






V 




0.0 


1000 


2.31109 


2.5 


1000 


18.7057 




642.233 


2.31117 




642.233 


18.7085 




316.228 


2.31119 




316.228 


18.7136 


0.5 


1000 


4.30196 


3.0 


1000 


23.9226 




642.233 


4.30221 




642.233 


23.9268 




316.228 


4.30248 




316.228 


23.9345 


1.0 


1000 


6.93473 


3.5 


1000 


29.7865 




642.233 


6.93529 




642.233 


29.7924 




316.228 


6.93609 




316.228 


29.8037 


1.2 


1000 


8.16826 


4.0 


1000 


36.2975 




642.233 


8.16900 




642.233 


36.3056 




316.228 


8.17011 




316.228 


36.3212 


1.5 


1000 


10.2122 


4.5 


1000 


43.4556 




642.233 


10.2133 




642.233 


43.4663 




316.228 


10.2150 




316.228 


43.4872 


2.0 


1000 


14.1356 


5.0 


1000 


51.2608 




642.233 


14.1375 




642.233 


51.2746 




316.228 


14.1406 




316.228 


51.3020 



TABLE 3. Condensate in Various Gauges for F(p) = 1 using WGTI 



Chapter 4 



Landau-Khalatnikov- 
Pradkin Transformations 
and the Fermion 
Propagator 

4.1 Introduction 

Looking for the gauge independence of the physical observables, we found that 
the Ward-Green- Takahashi Identity (WGTI) is a necessary condition for it, but 
not sufficient. We see, therefore, the need for the incorporation of other gauge 
invariance constraints for this purpose. We incorporate the Landa-Khalatnikov- 
Fradnin (LKF) transformations, [H H3 001 E3 EH] , which tell us the manner 
in which Green's functions change under a variation of gauge. Rules governing 
these transformations are far from simple. In counterpart, WGTI |341 1351 
are more simple, and therefore, they have been widely implemented. We can 
make larger this set of identities by transforming also the gauge parameter 
£ |59l lBT)| , and arrive to the Nielsen Identities (NI) . One advantage of these 
identities over the conventional Ward identities is that d/d^ becomes part of 
the new relations involving Green's functions. This fact was exploited in |611 
l(:>2| to prove the gauge independenc e of some quantities related to two-point 
Green's functions at the one-loop level and to all orders in Perturbation Theory, 
respectively. Since it is a difficult task to establish tha gauge independence of 
the physical observables in the Schwinger-Dyson Equations (SDE) study, NI can 
play a significant role in addressing this issue, along with the WGTI and the 
LKF transformations. However, in this chapter, we have focused only in the 
later. 

The LKF transformation for the three-point vertex is complicated and ham- 
pers direct extraction of analytical restrictions on its structure. Burden and 
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Roberts, [7j, carried out a numerical analysis to compare the self-consistency of 
various ansatze for the vertex, |19U37llS5| . by means of its LKF transformation. 
In addition to these numerical constraints, indirect analytical insight can be 
obtained on the nonperturbative structure of the vertex by demanding correct 
gauge covariance properties of the fermion propagator. In the context of gauge 
technique, examples are |64U65ll6*6*| . Concerning the works based upon choosing 
a vertex ansatze, references |18l 1371 l4"21 1431 ISlfl I58| employ this idea 1 . However, 
all the work in the later category has been carried out for massless QED3 and 
QED4. The masslessness of the fermions implies that the fermion propagator 
can be written only in terms of one function, the so called wavefunction renor- 
malization, F(p). In order to apply the LKF transform, one needs to know 
a Green function at least in one particular gauge. This is a formidable task. 
However, one can rely on approximations based on perturbation theory. It is 
customary to take F(jp) = 1 in the Landau gauge, an approximation justified 
by one loop calculation of the massless fermion propagator in arbitrary dimen- 
sions, see for example, |67| . The LKF transformation then implies a power law 
for F{p) in QED4 and a simple trigonometric function in QED3. To improve 
upon these results, one can take two paths: (i) incorporate t he information 
contained in higher orders of perturbation theory and (ii) study the massive 
theory. As pointed out in |2H1, in QED4, the power law structure of the wave- 
function renormalization remains intact by increasing order of approximation 
in perturbation theory although the exponent of course gets contribution from 
next to leading logarithms and so on. In |58j . constraint was obtained on the 
3-point vertex by considering a power law where the exponent of this power law 
was not restricted only to the one loop fermion propagator. In QED3, the two 
loop fermion propagator was evaluated in j'22l 1571 16*%| , where it was explicitly 
shown that the the approximation F(p) = 1 is only valid upto one loop, thus 
violating the transversality condition advocated in |56) . The result found there 
was used in |23| to find the improved LKF transform. 

In the present Chapter, we calculate the LKF transformed fermion propa- 
gator in massive QED3 and QED4 2 . We start with the simplest input which 
corresponds to the lowest order of perturbation theory, i.e., S(p) = l/ip'—m in 
the Landau gauge. On LKF transforming, we find the fermion propagator in an 
arbitrary covariant gauge. In the case of QED3, we obtain the result in terms 
of basic functions of momenta. In QED4, the final expression is in the form of 
hypergeometric functions. Coupling a enters as parameter of this transcenden- 
tal function. A comparison with perturbation theory needs the expansion of the 
hypergeometric function in terms of its parameters. We use the technique deve- 
loped by Moch et. al., for the said expansion. We compare our results with 
the one loop expansion of the fermion propagator in QED4 and QED3, |70U71| . 
and find perfect agreement upto terms independent of the gauge parameter at 
one loop, a difference permitted by the structure of the LKF transformations. 
We believe that the incorporation of LKF transformations, along with WGTI, 

1 A criticism of the vertex construction in lo(il was raised in J57| . 

2 In the context of gauge technique, gauge covariance of the spectral functions in QED was 
studied in 1531 1551 1551 . 
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in the SDE can play a key role in addressing the problems of gauge invariancc. 
For example, in the study of the SDE of the fermion propagator, only those 
assumptions should be permissible which keep intact the correct behaviour of 
the Green functions under the LKF transformations, in addition to ensuring 
that the WGTI is satisfied. It makes it vital to explore how two and three point 
Green functions transform in a gauge covariant fashion. In this Chapter, we 
consider only a two point function, namely, the fermion propagator. 



4.2 Fermion Propagator and the LKF Transfor- 
mation 

We start by expanding out the fermion propagator, in momentum and coordi- 
nate spaces respectively, in its most general form as follows : 

S,(p;f) = A{p;t)+x—j-= .j_ MM , (4.1) 

S F (x;0 = fX{x;Z)+Y{x;£), (4.2) 

where we explicitly write the gauge in which we specify the Fermion Propagator. 
Motivated from the lowest order perturbation theory, we take 

F(p; 0) = 1 and M(p; 0) = m . (4.3) 

Perturbation theory also reveals that this result continues to hold true to one 
loop order for the wavefunction renormalization. Eqs. (|4.1I4.2J) are related to 
each other through the following Fourier transforms : 

Sf(p;0 = / ' d*x(p-"S F {x;t) (4.4) 



Sf(x;0 = J J^Le-^'SFfaS) , (4.5) 

where d is the dimension of space-time. Let us recall that the LKF transforma- 
tion relating the coordinate space fermion propagator in Landau gauge to the 
one in an arbitrary covariant gauge reads l|2.98[) : 

S F (x;0 = Sf(z;0)e- l[Ad(0) - A ^ )] , 

and, also that H2.104[) 

/ r /d ri p.-ipx 

Taking ip to be the angle between x and p, we can write 
d d p = dpp d - x sin d ~ 2 ipdipCl d _ 2 , 
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where n d _ 2 = 2 n^-^/T ((d - l)/2). Hence 

A d (x) = -i£e 2 fi 4 - d f{d) / dpp d - 5 di>sm d - 2 ije- ipxcos,p , (4.6) 



where /(<f) = f2<i-2/(27r) d . Performing angular and radial integrations, we arrive 
at the following equation 

With these tools at hand, the procedure now is as follows : 

• Start with the lowest order fermion propagator and Fourier transform it 
to coordinate space. 

• Apply the LKF transformation law. 

• Fourier transform the result back to momentum space. 



4.3 Three Dimensional Case 

Employing Eqs. (|4.1I4.2I4.3I4.5[) . the lowest order three dimensional fermion 
propagator in Landau gauge in the position space is given by 

, , e- ,nx (l + mx) , . 

X( X ;0) = (4 .8) 

Y{ X ;0) = 



4irx 

Once in the coordinate space, we can apply the LKF transformation law using 
expression (|4.7|l explicitly in three dimensions : 

A 3 (*) = , (4.9) 

where a = e 2 /4ir. The fermion propagator in an arbitrary gauge is then 

S F {x; = S F (x; 0)e- (Q5/2)x . (4.10) 
For Fourier transforming back to momentum space, we use 
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Performing the angular integration, we get 



MP;0 = / dx sinpxe- {m+a ^ l2)x , (4.12) 

P Jo 

B(p;0 = - / -4{l + mx) [pxcospx-smpx}e~ im+a ^ 2)x ,(4.13) 



and the radial integration then yields 

4m 

4p 2 + (2m + a£) 



^(p;0 = . ,,2 ( 4 - 14 ) 



fl feO = - 4 f 2 4 - " + ffi + |i arctan [2p/(2m + afl] ■ (4.15) 
Ap^ + (2m + at,) z 2p 

We can now arrive at the following expressions for the wavefunction renormali- 
zation and the mass function, respectively : 

F(p;£) = arctan [2p/(2m + a£)] 1 2 P( 4 P 2 + a2 ^ 



2p L ^ /v S/J 0(p;£) 

a£(4p 2 +a£(2m + a£)) arctan [2p/(2m+a£)] 
^0 



(4.16) 



Mfofl = J^, (4.17) 

(pip; 

where 

0(p;£) = 2p(4p 2 + a£,(2m + a^))~a^(Ap 2 + (2m+a£,) 2 ) arctan [2p/(2m + a^)] . 

(4.18) 

In the massless limit, we immediately recuperate the well-known results : 

Fnm (p; = 1 - ^ arctan ^ , (4.19) 
M nm { P ;0 = 0. 

In the weak coupling, we can expand out Eas. l|4.16l4.17|l in powers of a. To 
0(a), we find 



M (p; £) = m 



F{p;0 = 1H q K™ 2 — P 2 ) arctan [p/ml — mp] , (4.20) 



2p : 



1 + —j {(to 2 + p 2 ) arctan [p/m] — mp} 



2p 3 

Let us compare these results with the ones obtained in |71j 

2p 



(4.21) 



Fr-ioopCp; = 1 + ^4 [(™ 2 ~ P 2 ) arctan [^/ m ] ~ TO ^] ' ( 4 - 22 ) 



-Mi_i oop (p;£) 



1 + ^-j { [£(to 2 +p 2 ) +4p 2 ]arctan[p/m] -£mp} 



(4.23) 
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We of course only expect the results to be in agreement upto a term propor- 
tional to a, as allowed by the structure of the LKF transformations. There is 
no such term in Eq. I]4.22[l . Therefore, the agreement is exact. Eq. I|4.21[) and 
Eq. (|4.23(l become identical only after we subract out the non- vanishing term in 
the Landau gauge from Eq. 1)4. 23J) to write out the subtracted mass function at 
one loop as : 

■Mf-ioopfeO = m 

One can numerically check that without the above mentioned subtraction, Eqs. H4.21I4.23JI 
approach the same value only in the large momentum regime. 



1 + — 4r {(to 2 + p 2 )arctan[p/m] — rap] 
2p 6 



(4.24) 



4.4 Four Dimensional Case 

Employing Eqs. H4.ll4.2l4.3l4.5fl . the lowest order four dimensional fermion pro- 
pagator in coordinate space is given by 

o 

TO 

X ^°) = -^2^2(tox), (4.25) 

9 

TTl 

Y(x;0) = --—K x {mx) : (4.26) 

47T Z X 

where K\ and Ki are Bessel functions of the second kind. In order to apply the 
LKF transformation in four dimensions, we expand Eq. (|4.7|l around d = 4 — e 
and use the following identities 



x e = l + elna; + 0(e 2 ) , 

to obtain 

(4.27) 



A 4 (x) ■ i( " 



2 

- + 7 + 2aifix + 0(e) 



167T 2 - e / 2 

Note that we cannot write a similar expression for A4 (0) because of the presence 
of the term proportional to In a;. Therefore, we introduce a cut-off scale x m i n . 
Now 

/ 2 \ v 

Ai{x m in)~ Ai(x) = -iln(-^— ) , (4.28) 

V X min J 

where v = a£/47r. Hence 

S F (x;Z) = S F (x;0)(^—) . (4.29) 

V X min / 

For Fourier transforming back to momentum space we use the following expres- 
sions 
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On carrying out angular integration, we obtain 



= / d A xp-xe tp - x X(x;0 



A(p;0 = x. 



2;/ 



/• oo 

-m 2 / <ixx' 
Jo 



dxx 



-2v+l 



K\ (mi) Ji (pre) 



if 2 (rec) J2 (px) 



The radial integration then yields 



A( P ;0 



1 1 m 
A 2 " 



2\ " 



r(l-i/)r(2-i/) 2 Fi l-i/,2-i/;2; 



where we have identified 2/a 



A. The above equations imply 



(4.31) 
(4.32) 

(4.33) 



B{P;0 = -^(tf) r(l-^)r(3-^) ^(1-^,3-^3;-^) .(4.3J) 



r(i-iz) 



2m 2 T(3-z/) 2 Fi(l-i/,3-i/;3;-p 2 /m 2 ) V A 2 



2\ " 



4m 2 r 2 (2 -i/) a-F? I 1 - 2 - 1/; 2 



+ P 2 T 2 (3-i/) 2j F\ 2 I l-i/, 3-1/; 3; 



M(p;0 = 



2m 2 Fi (1-v,2-l>; 2; -p 2 /m 2 ) 
(2-iz) 2 Fi(l-;/,3-i/;3;-p 2 /m 2 ) 



(4.35) 



(4.36) 



Eqs. I|4.35I4.36|I constitute the LKF transformation of Eqs. (|4.3|) . We shall now 
see that although Eqs. 1)4. 3JI correspond to the lowest order propagator, their 
LKF transformation, Eqs. (|4. 3514. 36(1 . is nonperturbative in nature and contains 
information of higher orders. 

4.4.1 Case a = 

Let us switch off the coupling and put a = which implies v — 0. Now using 
the identity 

2 F 1 (l,2;2;-p 2 /m 2 ) = 2 Fi(l, 3; 3; ~p 2 /m 2 ) = (1 + p 2 /™ 2 )" 1 , (4.37) 
it is easy to see that 

F(p; £) = 1 and M(p;£) = m , (4.38) 
which coincides with the lowest order perturbative result as expected. 
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4.4.2 Case m » p 

In the limit m >> p, the hypergeometric functions in Eqs. (|4.35I4.36|I can be 
easily expanded in powers of p 2 /m 2 , using the identity 



2-Fi a, P\ 7; o 



a/3 p 2 / p 2 



7 rrr 



+ 



(4.39) 



Retaining only 0(p 2 /m 2 ) terms, we arrive at : 

F(p;0 



r(i-z/)r(2-i/) z^to 2 

A 2 



(1-1//2) 
1+^ 



1- 



5zA p 



8 m 2 



-O 



M(p;0 = 



(1-I//2) 



(4.40) 



Now carrying out an expansion in a and substituting = a£/4-7r, we get the 
following 0(a) expressions : 



F( P ;0 



1 



47T 

m il + 



2 7 - 
8vr 



2p 2 
3m 2 



In 



A 2 



3m 2 



(4.41) 
(4.42) 



Let us now compare these expressions against the one-loop perturbative evalu- 
ation of the massive fermion propagator, see e.g., |70j : 



i*i-ioo P (p;£) 
Mi-i oop (p;0 
where 



47T 

am 

7T 



o< (1 -r 



£ m 2 
4 p ; 



2 



(4.43) 
(4.44) 



L = 
C = 



In 1 



7 — In 7r — In 



Knowing the fermion propagator even in one particular gauge is a prohibitively 
difficult task. Therefore, Eqs. (|4.3|l have to be viewed only as an approximation. 
For the wavefunction renormalization F(p;0), this approximation is valid upto 
one loop order, whereas, for the mass function, it is true only to the lowest 
order. Therefore we cannot expect the LKF transform of Eqs. I|4.3|l to yield 
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correctly each term in the perturbative expansion of the fermion propagator. 
However, it should correctly reproduce all those terms at every order of expan- 
sion which vanish in the Landau gauge at 0(a) and beyond. Therefore, we 
expect Eq. (|4.43|1 to be exactly reproduced and Eq. (|4.44|l to be reproduced 
upto the terms which vanish in the Landau gauge at O(a). After subtracting 
these terms, the resulting subtracted mass function is : 



Mf_ loop ( P ;0 



47T 



1+^(1-*) 



(4.45) 



In the limit m — » oo, the wavefunction renormalization acquires the form 



47T 



2p 2 



3m 2 



while the subtracted mass function is 



8n 



1 



3m 2 



(4.46) 



(4.47) 



The last two expressions are in perfect agreement with Eqs. (I4.41I4.42[1 after we 
make the identification : 



C7i £ -> 2 7 + In 



A 2 



(4.48) 



4.4.3 Case of Weak Coupling 

The case m » p is relatively easier to handle as we merely have to expand 
2-Fi(/3, 7; S] x) in powers of x and retain only the leading terms. If we want to 
obtain a series in powers of the coupling alone, we need the expansion of the 
hypergeometric functions in terms of its parameters (3 and 7. We follow the 
technique developed in One of the mathematical objects we shall use for 

such an expansion are the Z-snvas defined as : 

Z(n;m 1 ,...,m k ]x 1 ,...,x k )= 22 • ■ ■ ^ ■ ( 4 - 49 ) 

n>i 1 >i 2 >...>i k >0 1 k 

For xi = . . . = Xk = 1 the definition reduces to the Euler-Zagier sums, [73173] : 

Z(n;mi, . . . ,m k ; 1, . . . , 1) = Z mu _ <mh {n) . (4.50) 

Euler-Zagier sums can be used in the expansion of Gamma functions. For 
positive integers n we have |69| : 

r(n + e) = T(l + e)r(n) [l + eZi(n - 1) + . . . + e n_1 Zn.„i(n - 1)] . (4.51) 

The first sum Z\(n — 1), e.g., is just the (n — l)-th harmonic number, H n _i, of 
order 1 : 

71 — 1 ^ 

Z 1 (n-l) = J2- = H n - 1 . (4.52) 
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With these definitions in hand, we proceed to expand a hypergeometric function, 
2-Fi(l + £j 2 + e; 2; a;), as an example, assuming |x| < 1 : 



2 Fi(l + e,2 + e;2;2;) 



1 



T(2) 



E 



n=l 
= 1 + 



r(l + e)r(2 + e) 

r(l + e + n)r(2 + £ + n) x r ' 
r(2 + n) ~n\ 



E 



(l + e)r 2 (l+e) 

(1 + £ + n)(e + n) 2 T 2 (e + n) x n 
r(2 + n) nT 



(4.53) 



(4.54) 



Employing Eq. Ij4.51|l . we can expand the last expression in powers of e to any 
desired order of approximation. We shall be interested only in terms upto O(a). 



2 F 1 (l + e,2 + e;2;x) 



oo oo oo o 

Z — ' ^ ^n(n+l) 

n—l y ' 



n=l 



-2eJ2H n 



-ix 



(4.55) 



Performing the summations, we obtain 

1 



2 i ;l i(l+£,2 + £;2;a;) 



l-x 



I : { 1+ ln(l - 



Similarly, 



(4.56) 



2 Fi(l + £,3 + £;3;x) 



l-x 



1 3 1 
-e<- + -- 



x 2 1 — x 



1+x 



l-x 



In (1 - x) \ . (4.57) 



Substituting back into Eqs. (|4.35I4.36|I and identifying e = —v, we obtain 



F(p;0 



4tt 

ill 



27 - In 



A 2 



(1-L) 



(4.58) 



1 + 



which matches exactly onto the one loop result of Eqs. (|4. 4314. 4511 after the 
same identification as before, i.e., (|4.48|) . Therefore, we have seen that the LKF 
transformation of the bare propagator contains important information of higher 
orders in perturbation theory. 

LKF transformations, being non perturbative in nature, tell us the non per- 
turbative way in which the Fermion Propagator |74j and the Vertex transform 
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under a gauge variation. Most of the works on SDE violate these transfor- 
mations, with the exception of few of them which include the corresponding 
transformation to the Fermion Propagator at one-loop level. Since in Pertur- 
bation Theory all of the gauge invariance constraints are valid order by order, 
we can exploit this fact to construct a vertex ansatz such that in its pertur- 
bative expansion, it automatically satisfies, along with its associated Fermion 
Propagator, their corresponding LKF transformation, as we will do in the next 
Chapter. 



Apendix 



Most of the integrals involved in this Chapter are listed below for a quick refe- 
rence [SS1|7S] : 



£ <# sin d - 2 V cos ip e - tpxcos ^ = -i^r^y^T^' 1 



2 l-d/2 a d/2- 

For the three dimensional case, the needed integrals are : 

2 sinpx 



Jd(px) ,(4.59) 



(4.60) 



dJB sm9 e - ipxcosf> 



dO cos 9 sin 6 e 



— ipxcosU 



= 2i 



px 

cos px sin px 



dp — 5^ 

yr + mr) 



cospx smpx 
px {px) 2 

p sin px 



px (px) 2 
tt (1 + mx) 



dp- 



lo (p 2 + m2 ) 

— f ^e~ ax \px cospx — sin px] 
PJo x 

1 f°° dx 

- — e ax \px cos px — sin px] 
V Jo x 

dxsmpxe' {m+a ^ 2)x 



2 6 ~ 



a p 
— 1 H — arctan — , 
p a 



(4.61) 
(4.62) 
(4.63) 
(4.64) 
(4.65) 



arctan — , (4.66) 

a 



a 2 + p 2 p 



(m + a£/2) +p 2 



(4.67) 



For the four dimensional case, we used the following integrals in particular : 

[ ddsm 2 9e- ipxcose = — JApx) , (4.68) 
Jo PX 
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L dp WTm^ = ¥r^TT) K "-» {mx) > (469) 

X dxx = 2Wr( i +l/) r ( 2 J 

r (4.70) 



x i F 1 [ , 



Some of the series used in our calculation are as follows : 

f>-i*» = (4.71) 

y " + 1 r » 2 + .T (l + :r 2 )m(l-x) 

^-n(n + 2) ~ 4a; 2x 2 ' 1 J 

Y I x n 2-£ (l- g )]n(l-») (4?3) 



^ (n+l)(n + 2) 



Chapter 5 

Constructing the Vertex 



5.1 Introduction 

In studies on the Dynamical Breaking of Chiral Symmetry (DBCS) in QED3 
in the quenched and unquenched approximations [7J E] El El IHH1 EH1 UJ3 
C3 ESI E3 EDI ED IH2- we look for gauge independent physical observables. 
For that purpose, we have seen previously the necessity of including the Ward- 
Green- Takahashi Identity (WGTI) as well as the Landau-Khalatnikov-Fradkin 
(LKF) transformations. We know that, in Perturbation Theory, these gauge 
invariance constraints are satisfied order by order, thus, the vertex should be 
modified in every order of approximation. This fact has been exploited by, for 
example, [22U231 1571 168j . In this Chapter we carry out the constructionn of the 
non perturbative Vertex based on its perturbative counterpart, following the 
work . 

Making use of the WGTI, that relates the Vertex to the Fermion Propagator, 
one part of the Vertex, called longitudinal, can be expressed in terms of the said 
propagator 19 . We perform the evaluation of the Fermion Propagator at the 
one-loop level and hence we determine the longitudinal vertex at the same order. 
We also calculate the Full vertex at one-loop, and a mere substraction of the 
longitudinal part yields the transverse part, which is not fixed by the WGTI. 
According to the choice of Ball and Chiu, later modified by Kizilersii et. al. 
|7U|. the transverse vertex can be expressed in terms of 8 independent spin 
structures. Vertex should be free from kinematical singularities. Ball and Chiu 
choose the basis in such a way that the coefficient of every element of this basis 
is independently free from kinematical singularities in Feynman gauge. It was 
later shown by Kizilersii et. al. [21] that a similar calculation to the one of Ball 
and Chiu in an arbitrary covariant gauge does not have the same nice feature. 
Consequently, they proposed a modified basis whose coefficients are free from 
kinematical singularities in an arbitrary covariant gauge. The calculation in the 
present Chapter confirms that all the vectors of the modified basis also retain 
this feature for massive QED3. The final result for the transverse vertex is 
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written in terms of basic functions of the momenta in a form suitable for its 
extension to the non-perturbative domain, following the ideas of Curtis and 
Pennington pi?] . 

Using pcrturbative constraints as a guide, we carry out a construction of the 
non-perturbative vertex, which has no explicit dependence on the coupling a. 
This vertex has an explicit dependence on the gauge parameter £. For practical 
purposes of the numerical study of DCSB, we also construct an effective vertex 
which shifts the angular dependence from the unknown Fermion Propagator 
functions to the known basic functions, without changing its perturbative pro- 
perties at the one-loop level. We believe that this vertex should lead to a more 
realistic study of the dynamically generated masses through the corresponding 
SDEs. 



5.2 Longitudinal and Transverse Vertex to One- 
Loop 

5.2.1 The Fermion Propagator 

One-loop Fermion Propagator can be obtained by evaluating the graph in Dia- 
gram (7). 



-l 

• 

p 




Diagram 7 : One-loop Correction to the Fermion Propagator. 
This graph corresponds to the following equation : 

/rl 3 k 
p/^WfAj,(,), (5.1 

where q = k — p. The bare Fermion and Photon Propagators are, respectively : 

1 



s°Ap) = 



f > — m 

AJU?) = -feV + K-l)?*]/^ (5-2) 



where m is the bare mass of the fermion and Sp(p) represents the full propaga- 
tor, defined in eq. (|3.3|l . 

Taking the trace of Eq. 1)5. 1|) . having multiplied it with p 1 and with 1 re- 
spectively, one can obtain two independent equations. On simplifying, these 
equations can be written as : 



m - 1+ '^fwkw^) l(k2+p2}k - r - 2kV] ,5 - 3) 
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M( P ) 
F(p) 



/d 3 k 



) 3 q 2 (k 2 — to 2 ) 



(5.4) 



On Wick rotating to the Euclidean space and carrying out angular and radial 
integrations, we arrive at : 



1 



F(p) 
M{p) 
F(p) 



= [m-(TO 2 +p 2 )/(p)] , 

= m [1 + <*(£ + 2) J(p)] , 



(5.5) 



where we have used the simplifying notation I(p 2 ) = (1/y — p 2 ) arctan —p 2 /m 
Equations ((3.3(1 and ((5.5(1 form the complete Fermion Propagator at one loop. 



5.2.2 Longitudinal Vertex to One Loop 

We take the longitudinal part of the Vertex as the one that satisfies the WGTI, 
i. e., the Ball-Chiu Vertex, eq. 1(3.27(1 . On substituting eq. ((5.5(1 in the said 
expression, we obtain : 



r£(=r£ c ) = 



1 + T <ri 



(5.6) 



where 



01 
03 



k 2 
1 



I{k) 



k 2 — p 2 



k 2 



TO [I(k) - I(p)} 



TO 2 + p 2 



■7(A) 



I(p) — Tu- 



rn + p 2 



k 2 p 2 

I(p) + TO 



k 2 - P 2 ^ 



k 2 p 2 



(5.7) 



Eqs. ((5.6(1 and 1)5.7(1 give the longitudinal part of the fermion-photon vertex to 
one loop for the massive QED3. 

5.2.3 Full Vertex to One Loop 

Full Vertex can be obtained from Diagram (8) : 




Diagram 8 : One-loop Correction to the Vertex. 




so 
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and can be expressed as : 

r"(fc,p) = 7 M + A" . (5.8) 
Using the Feynman rules, A' 1 to 0(a) is simply given by : 

-ieA»= f 7 ^H e7 a )iSS.( P -«;)(-iey')<SS.(*-«')(-iey S )<A^(«;) ) 

(5.9) 

where the loop integration is to be performed in Minkowski space (in Euclidean 
space, these definitions are modified by the correspondig factors of i). A p can 
be expressed as : 



A M = - 



2 7T 2 



mV] 



^iVTa + 7"7 ,y 7 Al ^7a + 6m9""] + 7 Q 7 iy 7 M 7 A 7a^ I 



(2) 



+ (£-1) 



y^(0) _ [ 7 "jfy + 7^7" + 2m 3 H J« 



(5.10) 



+ [7" jf 7" # 7 A + m( 7^ i> 7 M 7 A + l v 7 P # 7 A ) + m 2 7^ 7^ 7 A ] 7^ 



where the integrals j( ), J^, J$?, 7<°), ijP and 7$ are : 



K<® = 




d 3 w 




1 






/ 

jm 


[(p — w) 2 — to 2 ] \{k — w) 2 - 


- m 2 




J(°) = 


j 


d 3 w 




1 






Jm 


w 2 [(p — w) 2 


— TO 2 ] [(k — w 


2 - 


TO 2 ] 


J (1) - 




d 3 w 




Wfj, 








j 

JM 


w 2 [{p — w) 2 


— to 2 ] [(k — w 


2 - 


TO 2 ] 


7(2) _ 




d 3 w 




W^Wu 






s 

JM 


w 2 [(p — w) 2 


— to 2 ] [(k — w 


2 - 


TO 2 ] 


7(0) = 


j 


d 3 w 




1 






JM 


w 4 [(p — w) 2 


— to 2 ] [(k — w 


2 - 


TO 2 ] 


7« = 


j 


d 3 w 




«V 






JM 


w 4 [{p — w) 2 


— to 2 ] [(k — w 


2 - 


TO 2 ] 




j 


d 3 w 










JM 


w 4 [{p — w) 2 


— to 2 ] [(k — w 


2 - 


TO 2 ] 



(5.11) 



We evaluate these integrals following the techniques developed in |19H57ll^5l 
I70| . The results are tabulated below, employing the notation A 2 = (k-p) 2 — k 2 p 2 
and X Q = (2/i7r 2 )A(°) for X = I,J,K. 
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d d 



w 



The Scalar Integral 

In arbitrary dimensions, 

j(o) = f 

J [(k — w) 2 — m 2 ][(p — w 2 ) — m 2 ]w 2 

Employing Fcynman parametrization, 

J(°) = T(3) f d d w f dai f da 2 f da 3 

J Jo Jo Jo 

S(a s -1) 



(5.12) 



with 



We define D as 



[ai[(fc — w) 2 — m 2 ] + a 2 [(p — w) 2 — m 2 ] + a3W 2 ] 3 
a s = cti + a 2 + a 3 . 



(5.13) 
(5.14) 



D 



ai[(k — w) 2 — to 2 ] + a 2 [(p — w) 2 — to 2 ] + a 3 w 2 
a s w 2 — 2(a\k ■ w + a 2 p ■ w) + ai(k 2 — to 2 ) + a 2 (p 2 — to 2 ) 
, 2{a 1 k + a 2 p)-w {a x k + a 2 p) 2 



a 



w 



+ 



at 



(aik + a 2 p) 2 ai(k 2 - to 2 ) + a 2 (p 2 - to 2 ) 



at 



= a. 



at 



+ 



a\k + a 2 p\ ' 



a., 



Let 
Then, 

J(0) 



a s / 

a s [ai(fc 2 - to 2 ) + a 2 (p 2 - to 2 )] - {a\k + a 2 p) 2 



a\k + a 2 p 
w ► w . 

a s 



2 / da\ / da 2 I da 3 
Jo Jo Jo 

I 



6{a s - 1) 



d d w 



(5.15) 
(5.16) 

(5.17) 



\w 2 +a s 2 {a s [ai(fc 2 — m 2 )+a 2 (p 2 — to 2 )] — (aifc+a 2 p) 2 }] 
Making use of the formulas 



3 ■ 



r d d w 
J 

r d d w 

J (w 2 - s) n 



i , ,„ . d r (n — I) d_ 



(5.18) 
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we can write 



j(0) = _ i7r # r ( 3 _ff| / d ai / da 2 / da 3 (5.19) 



o Jo 



S(a s -1) 



as~ 3 {—a s [ai(k 2 — to 2 ) + a 2 {p 2 — m 2 )] + («ifc + a 2 p) 2 } 3_ 2 
At this point, we remind the Cheng- Wo Theorem (S3] : 
Cheng- Wu Theorem 1 // 



„1 n / n \ 

■'O i=l V i=l / 



(5.20) 



i/ien we can write 



/=y n^y n****! i_ z) a <) ^ ' ( 5 - 21 ) 

where the set of a has been split into two nonempty sets 

{a'} = {a u ...,a k }, {a"} = {a k+1 , . . . , o^} . (5.22) 
Making use of this theorem, 

= _ Z7r f r f3--J / da 3 da 2 dai (5.23) 

x S( ai -1) 

c*s~ 3 {— a s [ai(k 2 — to 2 ) + ct2(p 2 — m 2 )] + (a.\k + a 2 p) 2 } 3_ 2 

After a trivial integration over a± , we obtain 

D = (1 + a 2 + a 3 ) d ~ 3 [m 2 (l + a 2 +a 3 )(l + a 2 ) - fc 2 (l + a 2 + a 3 - 1) 

-p 2 (l + a 2 + a 3 )«2 + a 2 p 2 + 2a 2 fc • p] 3 ~^ . (5.24) 

Using — 2k ■ p = q 2 — k 2 — p 2 , 

D = (l+a2+a 3 ) d ^ 3 [m 2 {l+ a 2+a 3 )(l+a 2 )-a 3 k 2 -a 2 a 3 p 2 -a2q 2 } 3 ^^ . (5.25) 

Therefore, 



J(°) = -i**r fs-fi Hda. r 



2 



o (l + a 2 + a 3 ) d - 



' -(5.26) 



[m 2 (l + a 2 + a 3 )(l + a 2 ) — ct 3 k 2 ~ a 2 a 3 p 2 — a 2 q 2 ] 3 
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In our three-dimensional case 



da 3 



da 2 
1 



[m 2 (l + a 2 + a 3 )(l + a 2 ) - a 3 k 2 - a 2 a 3 p 2 - a 2 q 2 ] 2 



ITT 



Jo 



(5.27) 



where 



Jo 



I da 3 
Jo Jo 



da 2 



1 



[m 2 (l + a 2 + a 3 )(l + a 2 ) - a 3 k 2 - a 2 a 3 p 2 - a 2 q 2 }2 
Being simpler the integration over a 3 first, we take 

D = m 2 (l + a 2 + a 3 )(l + a 2 ) - a 3 k 2 — a 2 a 3 p 2 — a 2 q 2 
= a 3 [m 2 (l + a 2 ) - k 2 - a 2 p 2 ] + m 2 (l + a 2 ) 2 - a 2 q 2 
= a 3 [a 2 (m 2 - p 2 ) + (to 2 - k 2 )] + to 2 (1 + a 2 f - a 2 q 2 . 

Therefore, 

da 3 



(5.28) 



(5.29) 



Jo 



poo pc 

= da 2 

Jo Jo 

- 2 / 

Jo 



[a 3 [a 2 (m 2 —p 2 ) + (m 2 — k 2 )]+m 2 (l+a 2 ) 2 — a 2 q 2 } 2 



oc 

2 / da 2 



[a 2 (m 2 — p 2 ) + (to 2 — k 2 )] 
1 



Jo 



{a 3 [a 2 (m 2 — p 2 ) + (m 2 — fc 2 )] + to 2 (1 + a 2 ) 2 — a 2 q 2 }2 
1 1 



CK3=0 



da 2 



\ot 2 (to 2 - p 2 ) + (TO 2 - fc 2 )] [ m 2(l + a2 )2 _ a25 2]i 



(to 2 -p 2 ) Jo 



da 2 



1 m 2 — k 2 
a 2 + m2 _„2 



Let 



a 2 + 



to 2 - fc 2 
to 2 — p 2 



[to 2 (1 + a 2 ) 2 — a 2 q 2 ] 2 



da 2 = dz . 



Integration limits are transformed in the following way : 

TO 2 -fc 2 



a 2 = => z = 



TO 2 — p 2 



a 2 — > 00 z — > 00 . 



(5.30) 



(5.31) 



(5.32) 



Now we rewrite the integrand as : 



2/i 1 \2 2 

to (1 + a 2 ) — a 2 q 



m z 1 + z - 



to 2 — fc 2 x 



TO 2 — ]3 2 



TO 



Z — 



TO 2 — ]3 2 
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= m z + 



(m 2 — p 2 ) 2 m 2 — p 2 



2 i „2 ' 



-zq* + q 



i 2 -k 2 



m 2 — p 2 



2m 2 (k 2 -p 2 ) - q 2 (m 2 -p 2 )~ 
m 2 — p 2 

m 2 {k 2 — p 2 ) 2 + q 2 (m 2 — k 2 )(m 2 — p 2 ) 



= m 2 z 2 + z 



+ 



(m 2 — p 2 ) 2 



= cz + bz + a = R , 



(5.33) 



where 



c = 
b = 

a = 



m 



2m 2 (k 2 - p 2 ) - q 2 (m 2 - p 2 ) 
m 2 — p 2 

m 2 (k 2 - p 2 ) 2 + q 2 (m 2 - k 2 )(m 2 - p 2 ) 
2 _ ^2\2 



X 



(m 2 — p 2 ) 



(m 2 — p 2 ) 2 



Consequently, 



Jo = - 



(5.34) 



dz 1 



1 2a+bz 
arctan ■ 



m 2 -p 2 7 m2- fc 2 z ^ (m 2 -p 2 )^E 2y^EVR 



We consider 

1. Evaluation at infinity 
2a + bz 



lim 



m 2 (k 2 - p 2 )[2m 2 - k 2 - p 2 ] + \ 
2my/—x( m2 ~ k 2 ) 



(5.35) 



(5.36) 



2. Evaluation at the lower limit 



2a + bz 



2 _ — j 



-X + m 2 (k 2 - p 2 ){2m 2 - k 2 - p 2 ) 
2m^/—x{m 2 — p 2 ) 



(5.37) 



and then, 



Jo = 



arctan 



m 2 (k 2 - p 2 )(2m 2 - k 2 - p 2 ) + \ 
2m^/^x(m 2 — k 2 ) 



— arctan 



-% + m 2 (fc 2 -p 2 )(2m 2 - k 2 - p 2 ) 
2m^x(m 2 - p 2 ) 



(5.38) 
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We take 



a = \(jp — w) 2 — to 2 ] 



b = [(k- w) 2 - m 2 ] . (5.48) 

We define now 



D = z[(p- wf - m 2 } + (1 - z)[(k - w) 2 - m 2 } 

= z(p 2 - 2p ■ w) + (1 - z)(k 2 ~2k-w) +w 2 -m 2 



Therefore, 



Using eqs. I|5.18[) . we arrive at 



= l\z{-l) 2 i^ J ^[m 2 -q 2 z{l-z)X 



o 

.2 / ,ir 2 2 



itt / dz[m - QT Z (l - z)] 2 



(I 



2in 2 -q 2 
: arctan ■ 



g 2 V 4m 2 



w 2 / I 



(5.49) 



We perform now the following change of variable 

w — > w — w — k(l — z) . (5.50) 

Then, we obtain 

D = zp 2 - 2zp ■ w' - 2p ■ kz(l -z) + k 2 z(l - z) + w' 2 - m 2 . (5.51) 
A second change of variable 

w — > w = w ~ zp , (5.52) 

implies 

D = w 2 + (k -p) 2 z(l - z) - m 2 

= w 2 +q 2 z(l- z) -to 2 . (5.53) 



K(0) = J! dz J ^WTq^hzT^Y 2 ■ (5 ' 54) 



which is our final expression. 
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The Tensor Integral 

We have 

jW = / d 3 w — . (5.56) 

^ J w 2 \{p — w) 2 — TO 2 ] [(k — w) 2 — TO 2 ] 

We write this integral in its most general form as : 

- 2 

Jjp = % ~Y {kMk,p)+p^J B (k,p)} . (5.57) 
Contracting with fc M and p^, we obtain the following system of equations : 

WJjP = ^[k 2 J A + (k-p)J B ] 

l?4 l) = l -^-[(k-p)JA+p 2 J B }. (5.58) 

On the other hand, 

(p 2 -TO 2 ) f. 1 



P 



W 2 [(p — w) 2 — TO 2 ] \{k — w) 2 — TO 



21 



2 J [(k — w) 2 — m 2 ]w 2 
If,, 1 



2 J [(p — w) 2 — to 2 ] [(A; — w) 2 — to 2 ] 
= ^^J(°) + \KV>- 5/(0,1,1). (5.59) 

1(0, 1, 1) will be dinned as a Master Integral when we calculate 1^ , and we use 
the relation : 

p • w = —(p 2 + w 2 — (p — w) 2 — to 2 + m 2 ) . (5.60) 
In a similar way, we also have that : 



k »jW = j(0) + 1^(0) _ l J(lj0j 1} . ( ,, n , 



We can now solve the system of equations (|5.58|l . On solving, we find 
JA(k,p) = 



-J2 j [P 2 (fc 2 -fc-p)-m 2 (p 2 -fc-p)] ^+k-pl(k) 

-p 2 Hp) + \ip 2 - k • p)/(g/2)| , 
J S (M = J^(p,fc). (5.62) 
Equations (|5.57|l and l|5.62|l form the final answer. 
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The Tensor Integral jffi 
We have 

J$ = J W 2 [(p _ w y Z ^2] p _ w )2 _ m 2] ' ( 5 ' 63 ) 

We write this integral in its most general form as : 

J jl 2 J = {^~ Ko+ (^^"SWy) Jc+ (p^K + k^Pu-g^ ^^ — ^ Jn 

Contracting with p**, 

- ^{y^o+((fe-p)fe,-p,y) Jc (5.65) 

+ (v 2 k u + (k -p)p v -p v — ^ P ^ JD+^p 2 Pi/-Pi/y^ Jgj • 
Performing the remaining contractions we obtain : 

(5.66) 

Other possible contractions are : 

k»P»J$ = ^[^K + \k\k-p)J c+ {kV+ { -^f-)jD 

+ lp 2 {k-p)J E \ (5.67) 



and 



{k ■ p)p v - k£\ J E \ , (5.68) 



jUM r(2) 



3 , 

wich, after a second contraction, yield 



WVJ® - {y^o + ^ 4 Jc + ^ 2 (fc • p)Jd + ((* • p) 2 - ^) J*} 

(5.69) 

and 

P^ 2) - ^{^i^o + ^ 2 (fc-p)Jc+(fcV + i(fc-p) 2 )^ 

+ ^ P 2 (k-p)J E \ . (5.70) 
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On the other hand, 

p^J = jd* 



p ■ w w v 

w 2 [(p — w) 2 — m?][(k — w) 2 — m 2 



+ 



(p 2 - m 2 ) 



d 3 i 



I 



d 3 w 



w v 



w 2 [(p — w) 2 — m 2 ] \(k — w) 2 — m 2 
w v 



[(p — w) 2 — m 2 ] [(k — w) 2 — m 2 ] 



~ / d A w 



w 2 [(k — w) 2 — m 2 ] 



{ ^^^ + \Mk,p)-\E v{ k). 



Similarly 



fc M J(2) = ( fc2 - ^ 



(5.71) 
(5.72) 



'v. v 2 

The new integrals needed are E l/ (p) and A u {k,p). Before preceding with the 



w 2 [(p — w) 2 — m 2 ] 



(2) 

calculation of J^J , let us pay attention to them. First, 

E v {p) = J d 3 w 
Using Feynam Parametrization, 
E v (p) = dx J 



21 ' 



d A w- 



Let 

We take now 

such that 

and, consequently, 



x[(p — w) 2 — m 2 ] + (1 — x)w 2 

D = x[(p - w f - m 2 } + (1 - x)w 2 
= w 2 + xp 2 — 2xp ■ w — xm 2 . 

w = w + xp , 
D = w' 2 + p 2 x(l — x) — m 2 x , 

w' v +p u x 



(5.73) 
(5.74) 

(5.75) 
(5.76) 
(5.77) 



MV) = J^dxjd'w' — 

= p v I dxx(-l)'- 
Jo 



+ p 2 x(l — x) — m 2 x] 2 
«r* F p (2 ^ [-P 2 x(l -x) + m 2 x]i- 







m 2 p v j 


f dx 




iir 2 p v 


to — 




E v {k) . 





dx x 2 [— p 2 (1 — x) + m 2 



(to 2 - p 2 ) 



arctan 



(5.78) 
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The other integral we need to calculate is 

A v (k,p) = / d 3 W- ^ ^TTT 21 ' ( 5J9 ) 

J l(p — w) z — m l \ \{k — wY — m z \ 
After Fcynman parametrization, we write : 



A v {k,p) = J dz J d 3 wj : 



[z[(p — w) 2 — m 2 ] + (1 — z)[(k — wY — m 2 ]} 2 

(5.80) 

We define 

D = z[(p - w) 2 - m 2 } + (1 - z)[(k - wf - m 2 } . (5.81) 
After the change of variable 

w' = w-zp-(l-z)k, (5.82) 

we have 

D = w 12 + q 2 z{\- z)-m 2 , (5.83) 

and, therefore, 

Mk, P ) = fdzf d*w <+ P f + k ^- z) . (5.84) 
J J [w' z + q A z(l - z) - m 2 Y 

On simplifying, 

A v {k,p) = Cdzr, 2 P; z + k ^- z ) 
Jo \Jm 2 - q 2 z(l - z) 

= iTT 2 (k+p),j I dz- 
Jo 



lo \Jm 2 — q 2 z(l - z) 



m2{k+p) »-,rct,n\^= 1 -K^\k + p) v . (5.85) 



V /Z^2 V 4m 2 2 

Let us turn back to the claculation of jffl . Second contraction yields : 

P^jjS = (P2 7 2> TP'P) J a +P 2 Jb] + \v v A v -^E v (k) 
k'YJjii = [I ' T/" , ! —[l- 2 J A +{k-p)JB] + \k v A v -^-E v {k) 



(p 2 


-m 2 ) 


iir 3 


(p 2 


2 

— TO 2 ) 


2 

iir 3 


(k 2 


2 

— TO 2 ) 


2 

in 3 


(k 2 


2 

— TO 2 ) 


2 

in 3 




2 


2 



= V " 9 '" >- ^[k 2 J A + (k-p)JB} + lk»A„-^E v (p) 

P^jff = V " ~ m ' ^[(k-p)J A +p 2 J B } + ^A u -^EAp) -(5.86) 
We arrive then to the following system of equations : 

J c + ^2 (fc . f))jD+ 2 p4jB = a 
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2 -k\j c + h 2 (k-p)J D +(^k-p) 2 -^f^J E = 



h 2 {k-p)J c +[k 2 p 2 



{k-pf 



Jd + ^p 2 (k-p)J E 



= c 



(5.87) 



where 



ITT 



b = 



(5.88) 



Solutions to this system are 
Jc(k,p) 



J D (k,p) 



i^j [p 2 (k ■ p-2k 2 )-m 2 (k • p-2p 2 )} ^-p 2 (p 2 m 2 )^ 

+ k -jl (- 2 - fc2 ) m + \ ( k -P + P 2 ) ^9/2) - | , 

^|[fc 2 (3fc. p -p 2 )-m 2 (3fc. p -fc 2 )] ^ 
+ [p 2 (3fc-p-fc 2 )-m 2 (3fc-p-p 2 )] ^ 
-(m 2 - fc 2 )/(fc) - (m 2 -p 2 )I{p) -\{k+ p) 2 I(q/2) + 2m| , 



J E (k,p) = Jc(p,k) , 
which completes or calculation of . 

The Scalar Integral 

We have 

7 (0) = / d 3 w - 

J w A [{p — w) 2 — m 2 ][(k — w) 2 — m 2 ] 

In general, let us define the family of integrals : 



(5.89) 



^1,^2,^3) 



d d i 



' [(p - w) 2 - m 2 ] Vl [(k - w) 2 - m 2 }" 2 [w 2 ] V3 
We use the following Integration by Parts (IBP) identity : 



(g» ~ w)jj, 

[(p - w) 2 - m 2 ] Ul [(k - w) 2 - m 2 }^ [ w 2 ] 



- . 



(5.90) 



(5.91) 



(5.92) 
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After the differentiation we obtain : 



J d d w{ qi - m) 

I 



2v 1 (p-wY 2v 2 {k-wY _ 2v 2 w» c „ 3+1 



+ I d d w _J_±_ = o , 



A" 1 B U2 C v z 



since 



d 
dw. 



We have defined 

A= {p-wf -m 2 , B = (k-w) 2 -m 2 , C = w 2 

Therefore, 



(5.93) 



(5.94) 



(5.95) 



dl(i/i, 1/2,1/3) = Jd d w(q l - 



w 



2v 1 {p-w) ti 2i/ 2 (k- W y 2v 2 w^ 



+ 



A Vl+1 B U2 C V3 A" 1 B" 2+1 C V3 A Ul B U2 C U3+1 



. (5.96) 



Assuming that qi is defined such that q\ — p, q 2 — k and q3 = 0. Then, for 
i = 0: 



dl (1/1,1/2, 1/3) = / d d w(p-w) 



2u 1 (p-wY 2v 2 {k-wY 



2v 2 w^ 



Avi + l Qv 2 (Jv 3 Jl s y\'Qv- 2 -\-\(Jv 3 (Jv 3 + l 



. (5.97) 



Let us observe the following products : 
f. 

2(p — w) ■ (k — w) = (p— w) 2 + (k — w) 2 — (p — k) 2 — m 2 — m 2 + 2m 2 
= A + B - [{p- k) 2 - 2m 2 } . (5.98) 



2. 



2(p-w)-(-w) = (p-w) 2 -m 2 +w 2 -p 2 + m 2 = A + C-p 2 + m 2 . (5.99) 



3. 



(p - w) 2 = (p-w) 2 -m 2 +m 2 = A + m 2 . (5.100) 
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Substituting these products into eq. I|5.97[) . we obtain : 

2v 1 [A + m 2 } v 2 [A + B-[{p-kf 



cfl(i/i, 1/2,1/3) 



d d w 
v 3 [A + C- (p 



2m 2 



2 )] 



1v\ + z/ 2 + "3 



d d w 

A^B V2 C V3 
v 2 [{p — k) 2 — 2m 2 } 



A v \ B U2+1 C" 3 



2v\m 2 



v?, 



A" 1+1 B U2 C 3 
v 3 {p 2 - m 2 ) 



^iJJ^ + l^ j\l Jl ~ 1 B' y2 C U3 + 1 B"2 (7^3 + 1 

= (2i> 1 + v 2 + V3)I(vi,i>2 1 v 3 ,) + v 2 l{vi-\,v 2 + 2,v 3 ) 

+2m 2 v 1 I(v 1 + l,is 2 ,v 3 ) - v 2 [{p-k) 2 - 2m 2 ]/(iA L , ^2 + 1,^3) 
+^3/(^1-1,^2,^3 + 1) - v 3 {p 2 - m 2 )I{vi, v 2 ,v 3 + 1) (5.101) 

Rearranging the terms in such a way that those integrals with V\ + v 2 + v 3 
as the sum of their arguments are on the right hand side, while those with 
v\ + v 2 + V3 + 1 on the left hand side, we get : 

-2^im 2 /(^i + l, i/ 2 , V3) + v 2 [{p - k) 2 - 2m 2 \l(y ll t/ 2 + l, v 3 ) 

+v 3 (p 2 - m 2 )I(vi,v 2 , V3 + 1) = (2vi + v 2 + V3-d)I{vi,V2,v 3 ) 

+u 2 I{vi -l,u 2 + 1, v 3 ) + v 3 I{v x - 1, u 2 , v 3 + 1) . (5.102) 



For i = 2 we have 
dl{v 1 ,v 2l v 3 ) 



d d w(k — w)^ 

2v x (p-wY 2v 2 (k-wY 



2v 3 w^ 



A Vl+1 B" 2 C Vi A Ul B U2+1 C V3 A Ul B V2 C V3+1 
(Vl + 2v 2 + V3)I{vi,V2,V3) + v\I(y\ + l,v 2 - 1, V3) 
+v 3 I{vi,v 2 - 1,^3 + 1) 
-[{p - k) 2 - 2m 2 \v 1 l(v l + 1, i/ 2 , v 3 ) 
+2v2m 2 I{v\,V2 + l,v 3 ) 

^ 3 (fc 2 -m 2 )/K v 2 ,v 3 + l) (5.103) 



"i[(p- k ) 2 - 2m 2 ]I{v 1 + 1,^2,^3) - 2v 2 m 2 I{v l ,V2 + 1,V 3 ) 

+v 3 {k 2 - n{ 2 )l(y 1 ,V2, v 3 + 1) = (z/ 2 + 2i/ 2 + v 3 - d)I(y x ,V2,v 3 ) 



and for i = 3 : 
dl{vx,v 2 ,v 3 ) 



+ VlI(l>l + 1,V 2 - 1, l/ 3 ) + v 3 I{v Xl v 2l v 3 + 1) , 



d d w{-w^) 



(5.104) 
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2v 1 (p - wY 2v 2 {k - wY 2u 3 w^ 



or 



A" 1+1 B 1 " 2 C V:i A Vl B V2 + 1 C U3 A Ul B V2 C V3+1 

= {vi + V2 + 2v 3 )I(ui, v 2 , 2/3) + 2/1/(2/1 + 1, 2/2, v 3 - 1) 
+v 2 I{vi +V2 + I, v 3 - 1) 
-vi(p 2 -m 2 )/(^ 1 + 1, 2/2,2/3) 

-v 2 {k 2 -m 2 ) I \vx,v 2 + 1,1*) (5.105) 

v x {p 2 - m 2 )l(v x + l,v 2 ,v 3 ) + ^(fc 2 - ra 2 )/(z/i,!/ 2 + 1, 2/3) = 
(i/i + i/ 2 + 2z/ 3 - d)I(v±,v 2 , v 3 ) + vj{vi + 1, v 2 , v 3 - 1) 

+2/2/(2/1,2/2 + 1,2/3-!) ■ (5.106) 
Therefore, we need to solve the system of equations : 

- 2v 1 m 2 I{v 1 + 1, v 2 , 2/3) + v 2 (q 2 - 2m 2 )I(v 1 + v 2 + 1, 1/3) 

+K-s(p 2 - ™ 2 )^(^i, f2, 2/3 + 1) = a 
vi{q 2 - 2m 2 )I(y 1 + 1,2/2,2/3) - 2v 2 m 2 1{v x ,v 2 + l,v 3 ) 

+v 3 {k 2 -m 2 )I{v 1 ,v 2 ,v 3 + l) = b 
vi(p 2 - m 2 )I{v\_ + 1, u 2 , 1/3) + v 2 {k 2 - m. 2 )I(i/i, v 2 + 1, v 3 ) 

+0/(1/1,2/2,1/3 + 1) = c, (5.107) 

where 

a = (2i/ x + v 2 + 1/3 - d)i"(z/i + z/ 2 + 1/3) + 2/2/(2/1 - 1, J/2 + 1, 2/3) 

+i/ 3 /(z/i - l,z/ 2 , i/ 3 + 1) 
6 = (1/1 + 2z/ 2 + 1/3 - d)/(z/i + z/2 + 2/3) + ^1/(2/1 + 1, J/2 - 1, 2/3) 

+2/3/(2/1,2/2 - 1,2/3 + 1) 

C = (i/1 + f2 + 2l/ 3 - d)/(z/l + f2 + 2/3) + 2/1/(2/1 + 1, 2/2,2/3 - 1) 

+1/2/(2/1,2/2 + 1,2/3-1). (5.108) 
In obvious notation : 



a 



-2m 2 v Y v 2 (q 2 -2m 2 ) v 3 (p 2 - m 2 ) \ ( h 
Vl (q 2 -2m 2 ) -2m 2 v 2 v 3 {k 2 - to 2 ) J 2 j = I b 

v x {p 2 -m 2 ) v 2 {k 2 -m 2 ) / \ I 3 ) \c 



(5.109) 



which formally we can write as : 

SI = A . (5.110) 

Then, 

|S| = 2 I / 1 i/ 2 i/3[m 2 (fc 2 - m 2 ) 2 + m 2 (p 2 - to 2 ) 2 + (k 2 - m 2 )(p 2 - m 2 )(q 2 - 2m 2 )} , 

(5.111) 
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and besides 



\S\I 3 = ^ 2 [(q 2 -2m 2 ){k 2 -m 2 ) + 2m 2 (p 2 -m 2 )]a 
+viv 2 [(q 2 - 2m 2 )(p 2 - m 2 ) + 2m 2 (fc 2 - m 2 )]b 
+v 1 v 2 [4m 4 - (q 2 - 2m 2 ) 2 ]c . (5.112) 



For v 1 = v 2 = v- i = \, 



a = 7(1, 1,1) + 7(0,2,1) + 7(0, 1,2) 
b = 7(1, 1,1) + 7(2, 0,1) +7(1, 0,2) 
c = 7(1, 1,1) + 7(2, 1,0) + 7(1, 2,0) (5.113) 



with 7(°) = 7(1, 1, 2) and = 7(1, 1, 1). So, we need to calculate the integrals 

7(0, 2, 1) k = p 7(2, 0, 1), 7(0, 1, 2) k = p 7(1, 0, 2) and 7(2, 1, 0) k = P 7(1, 2, 0). For 
that purpose, we define the integral 

Ip n (k,p,m 1 ,m 2 ) = I d d wj- ^ gT^T? ^ 2i7T ■ (5.114) 

7 [(k — w) z - mi\P[(p — w) z - m^y 1 



Using Feynam parametrization, we know that 

0-*)" 't-i '—r. — (■ r »-ii-') 



1 _ r ("+P) f 1 dxx n-l (1 _ x)P -l L 

A n BP T(n)T(p)J ( ' [xA + (1 - x)B] n +P ' 



Let 



and 



^4 = [(p - w) 2 - m 2 ] , 

B = [(k - w) 2 - m 2 } , (5.116) 



7? = + (1 - x)B 

= x[p 2 + w 2 — 2p-w — m 2 ] + (1 — x) [k 2 + w 2 — 2k ■ w — m 2 ] 

= w 2 -2w ■ \px+k(l-x)]+p 2 x+k 2 (l-x)-m 2 ,x-rn 2 (l-x) . (5.117) 

We make the change of variable 

w' = w - \px + k(l - x)] , (5.118) 

such that 

w 2 - 2w ■ [px + fc(l - x)] = w' 2 - [px + k{l - x)] 2 . (5.119) 

In this way, D is rewritten as : 

D = w' 2 - [px + fc(l - x)] 2 +p 2 x + k 2 (l - x) - mix - m 2 (l - x) 

= w' 2 + p 2 x(l — x) + fc 2 x(l — x) — 2k- px(l — x) — m\x — m\{l — x) 
= w' 2 +q 2 x(l-x)-m 2 x-m 2 (l-x) . (5.120) 
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Therefore, 



T(n + p) f 1 



[ dxx n - 1 (l-x) p - 1 fd d w'— l — 
Jo J D n +P 



I pn (k,p,m u m 2 ) = r(n)r(p) 

x) 



T(n + P ) f 1 
T(n)T(p) 

i{ ir+V / 2 r(n + p-j/2) f\ 



Jo r ( n - 



d/2) 



P-lef-n-P 



with 



s = —q 2 x(l — x) + m\x + m\(l — x) . 
For the integrals of our interest, we help ourselves with 

d d w 



I 11 (k,p,m 1 ,m 2 ) = f J£p~Z 



w) 2 - m?,][(k - w) 2 - m\] 



m 2 T 2 



dxs 2 ' 



—in dxs 2 , 



where, in the last line, we set d = 3. So 

1 d 

7(2,1,0)= 7rzr -^ ::r hi(k,p,mi,m2) 



Now, 



and therefore, 



2m\ dm\ 
d 



ra\—ra2—ra 



dmi 



s 2 = — mi(l - x)s 2 



1(2 



dx- 



(1-x) 



[— q 2 x(l — x) + m 2 ] 2 



(5.121) 
(5.122) 



(5.123) 

(5.124) 
(5.125) 



m[4m 2 — q 2 



k = P 1(1, 2,0). 



On the other hand 

7(2,0,1) = 



_L d_ 

2m i dm\ 
2 r 1 



7ii(fc,p,mi,TO 2 ) 



k— mi— 0,m2— m 



2 i 



o [— p 2 x(l — x) + m 2 x] 2 

2 



m(m 2 — p 2 ) 



k = P /(0.2.1). 



(5.126) 



(5.127) 
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To solve 1(1, 0, 2), we take the integral 

r =: / d d w 



w 4 [(p — w) 2 — to 2 



(5.128) 



Then, 



■(p 2 - to 2 ) / d d w 



W 4 [(p — w) 2 — TO 2 



-- / d d l 



1 f d d i 



W 2 [(p — w) 2 — TO 2 ] 



= l -(p 2 - m 2 )I(l, 0,2) + |/(1,0,1) , 

since tha last term vanishes in the dimensional regularization scheme, 
to see that 



(5.129) 
It is easy 



1(1, 0,1) 



hi(k,p,m 1 ,m 2 ) 



k— mi— 0,m2— m 



in / \—p x(l — x) + TO x] 
Jo 



2iTT 2 



: arctan 



~P 



7(0,1,1) • 

To calculate explicitly we take in eq. (|5.115|l 

A=[(p- w) 2 - to 2 ] , B = w 2 , n = l, p=2 . 
So, in three dimensions, 



r 



r(3) 



r(i)r(2) y 



(1 - x) / d 3 ww; 1/ 



Let 
Then, 



w = w — px , 



1 



[Ax + B(l -x)} 3 



D 



Ax + B(l-x) 

w 2 + p 2 x — 2p ■ wx — xm 2 

12 i 2 2 2 2 

W — TO X 

It/ 2 +p 2 x(l — x) — TO 2 CC , 



and therefore 

r = 2 / dx(l -x) / (d 3 w)- 



p x 



w 2 +p 2 x(l — x) — m 2 x] 3 



(5.130) 
(5.131) 

(5.132) 
(5.133) 



(5.134) 



(5.135) 
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Using now H5.18fl . 



3 r(£) 3 

2p u I dxx(l-x)(-l) 3 int- F ^[-p 2 x(l-x)+m 2 x}2 



in 



r(3) 

p" / dxx~?(l - x)[-p 2 {\ - x) + to 2 



dx 



-P \ 2 P 



jx 2 - x* 1 [-p 2 (l - x) + 
2m 2 



to 2 — p 2 V to 2 



1 2 arctan ' 



So that 



7(1,0,2) 



p2 _ m 2 



Pu r - -i(i,o,i) 



2m 



(to 2 — p 2 ) 2 
7(0,1,2) . 



Finally, substituting into (|5.112|) 

/(°) = ~ { g 2 (m 2 + fc-p)j(°) + in 2 mL} 

where 

g 2 (fc 2 - m 2 ) - (fc 2 - p 2 )(fc 2 + to 2 ) 
(k 2 — to 2 ) 2 
q 2 (p 2 — m 2 ) + (fc 2 - p 2 )(p 2 + to 2 ) 



L = 



(p 2 — TO 2 ) 2 



This is our final expression. 

The Tensor Integral 1^ 

We have 

IP = I d 3 w 



I M w4 [iP ~ W ) 2 ~ m2 \ [(k — Vj) 2 — TO 2 ] 

We write this integral in its most general form as : 

_2 



4 X) = yMA(*.P)+Vsft?)] 



(5.136) 



(5.137) 
(5.138) 



(5.139) 



(5.140) 



(5.141) 

Contracting with fc M and p^ we obtain the following system of equations : 



k^ 



-[k 2 I A + (k-p)I E 



P^ 1} = — [(k-p)I A +p 2 I B ] 



(5.142) 
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On the other hand 



= 



(p 2 — TO 2 ) 



dh 



1 



; / d 3 i 



W A [(p — U}) 2 — TO 2 ] [(k — U}) 2 — TO 2 ] 
1 



[(k — w) 2 — m*\ir 



; / d 3 t 



21,„4 
1 



! [{p — Vj) 2 — TO 2 ] [(k — w) 2 — TO 2 



£^±W+\jV> -1/(0, 1,2), 



where we have used : 



1 



p ■ w = — (p 2 + w 2 — (p — w) 2 — m + m 2 ) . 



We also have that 

^ 2 " ' 2" 2 

We solve the system of equations (15.14211 , to find : 



fc/xjW = (fc 2 -m 2 ) j(01 + 1 J(0) _ 1 



J A (fc,p) = [fc-p(p 2 -TO 2 )- P 2 (fc 2 -TO 2 )] 



7(1,0,2) 



(5.143) 
(5.144) 
(5.145) 



mp 



mk ■ p 



(to 2 — p 2 ) 2 (to 2 — /c 2 ) 2 f 

I B {k,p) = I A (p,k) . 
Equations (|5.141f> and 15.146(1 form the complete solution. 

(2) 

The Tensor Integral I^J 

We have, 



(5.146) 



r(2) = 

1 \IV 



d 3 ~ 



l M w 4 [(p — w) 2 — to 2 ] [(k — w) 2 — m 2 ] 
We express this integral in its most general form as : 



(5.147) 



j(2) 



ITT 



9V 



k 2 \ ( 2(fc-p) 
Jq+ k^kv-g^— I c + PttK + k^Pv-g^ 



ID 



P 



+ \PnPv - J Ib 

Contracting with p^ we obtain : 

¥ LLV o 



(5.148) 



+ \(h ■ p)k v - pv— I c 



+ [p 2 k u + (k ■ p)p v -Pu 



2(k-p) 



lD+[p 2 Pu-p l ,hr)lE \ .(5.149) 
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Performing the remaining contractions, 



p v p^J - ^ {y J o + ((* • P f -^f)ic + \v 2 {k ■ p)I D + \ P H E } 

(5.150) 

and 

= ^l^Jo + (^\k.p))lc + (kV + ^-)ln 



3 

In the same way, 



l P 2 (k-p)I E ^ . (5.151) 



k"I$ = ^ jy Jo + ^ 2 fc^ c + ((fc • p)fc„ + fcV - |fc„(fc • p)) 7 D 

+ ((fc-p)?,-^)^} , (5.152) 
from where, after a second contraction, we obtain : 



k"k»I$ = | y Jo + Hfc% + ^ 2 (fc • p)I D + ((fc • p) 



2 fcV 



(5.153) 



and 

^3 



p-k^j = ^{^j Q+ i k \k.p)ic+(kv+i(k-pfy D 



+ lp 2 (k-p l ; . (5.1.14) 



On the other hand, 



^ J w 4 [{p - wf - m 2 ][(k - w) 2 - m 2 ] 

(p 2 -m 2 ) /■ 3 «;„ 



2 7 ^[(p-u;) 2 -™ 2 ]^-^) 2 -™ 2 ] 



-II 



d 6 w- 



[(p — w) 2 — m 2 ] [(k — w) 2 — m 2 



Similarly, 



w 4 [(k - w) w - m 2 } 
^ + \W ~\E v {k) . (5.155) 

= iJ ^- 1 ll 1) + \EA P ) . (5.156) 
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From the second contraction we obtain : 



pW$ 



(p 2 


-to 2 ) 


in 3 


(p 2 


2 

-TO 2 ) 


2 

Z7T 3 


(fc 2 


2 

-TO 2 ) 


2 

in 3 


(fc 2 


2 

-TO 2 ) 


2 

17T 3 



[(fc.p)/ A +p2/ B ] + -p"4 1 )-^-£7 w (fc) 

[fc 2 / A + (k ■ P )I B ] + \k v jM - y E v (fc) 



2 p-#i+A] + /4 1L y^(p) • (5.157) 



In this way, we arrive at the following system of equations : 

4 ,,. 2 

3 



(fc-p) 2 -^? 



Ic + ip 2 (k-p)I D + ^p 4 I E 



= a 



'-k i I c + -k 2 {k-p)I D + (fe-p) 



\2 



-fc 2 (fc-p)/ c + feV + 



(fc-p) 1 



= c, (5.158) 



where 



a = 



c = 



2 ., p 2 



Z7T° 



k ^m+p^m-itAj 



(5.159) 



Solution to the system are : 
/c(fc,p) 



^p 2 Jo + [P 2 (k ■ p 2k 2 ) - TO 2 (fc • p 2p 2 )] f 

pV - m r-f+(k- P - 2 P 2 ) f p 2 J -f ^ /(*) 

mfc • p | 



/a 



Ib 



Ja 



k 2 (m 2 — k 2 ) 



I D (k,p) = 



2A 2 



2fc • p J + [fc 2 (3fc -p-p 2 ) - TO 2 (3fc • p - fc 2 )] y 



+ [p 2 (3fc • p - fc 2 ) - TO 2 (3fc • p - p 2 )] ^ + (3fc • p - fc 2 ) 



^4 

2 



+ (3fc'P-p 2 ) 4r + I(k)+I(p) 



TO 2 — fc 2 TO 2 — p 2 



I E {k,p) = Ic(p,k) . 
This completes the calculation of the vertex to 0(a). 



(5.160) 
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5.2.4 Transverse Vertex to One loop 

We can subtract from the full vertex, Eq. I|5.11|) . the longitudinal vertex, Eqs. (|5.6(l 
and H5.7J1 . and obtain the transverse vertex to O(o). Let us recall that the trans- 
verse vertex T^,(k,p) can be written in terms of 8 basis vectors as follows : 

8 
i=l 

with the basis : 



with 



T>' 

rpfJ, 
J 2 

rp/l 

J :i 

rpfJ, 

J 4 

rp/l 

1 5 

rplJ, 

1 6 



rri/l 

1 S 



\p»{k-q)-W(p-q)] 
[p"(fc-g)-fc"(p •?)](/«+!<) 

q 2 jf + _ *J* _ p"] - 2{k - pYk x p v a Xv 



qv<3 

- 7 ^(fc 2 -p 2 ) + {k + P Y4 
-kk 2 - P 2 )ir(k' + p')-k» 



- 1 l *k l 'p x a ljX + k li p l -p li ft 



P»] + {k+pYk x p"a x „ 



(5.161) 



After a lengthy but straightforward algebra, the coefficients can be identified. 
We prefer to write these out in the following form : 



n{k,p) 



^2aij(k,p)I(lj 



k 2 p 2 



(5.162) 



where l\ = ?7 2 x/4, l\ = I 2 = k 2 , l\ = p 2 and I 2 = q 2 : /4. Functions rji, n 2 

and % have been defined above. Similarly, the factors gt are — gi = mA 2 g2 = 
2mA 2 g 3 = 2A 2 <74 = g 5 = 2mA 2 g 6 = A 2 gt = mg$ = m/4A 2 . The coefficients 



in the one loop perturbative expansion of the r,, Eq. H5.162f> . are tabulated 



below : 

a n (k,p) 
a 12 (k,p) 

ai3(k,p) 

a u (k,p) 
a\b{k,p) 
a 16 (k,p) 

a2i(k,p) 



-(£ + 2) ?/l (m 2 + k-p) 
an(p, k) 

(k 2 + k-p) 



4(£ 



(k 2 -p 2 ) 



ai 3 (p, k) 

-2(e + 2) 





-T)l< 



a 2 

^-m 4 + { (k ■ pf - (k 2 + p 2 )(k ■ p) + k 2 p 2 }m 2 
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a-22(k,p) 



a 2i {k,p) 
a2b{k,p) 



azi{k,p) 



-«-{ {k .pf + kV} 



+- 



2X 



g 4 m 8 - g 2 { (jfe • pf + 2(k 2 + p 2 )k ■ p - 5k V } 



+ ^ 2 (fc 2 +p 2 )A 2 m 4 



a2n{k,p) = toA < A: • p + 



+ {2(fc 4 + p 4 + k 2 p 2 )(k ■ pf - 7k 2 p 2 {k 2 + p 2 ){k ■ p) 2 
+10fcV fc • p - k 4 p 4 (k 2 + p 2 )}m 2 

+ l -k 2 p 2 q 2 {(k 2 + p 2 )(k ■ p) 2 - Ak 2 p 2 k ■ p + k 2 p 2 (k 2 + p 2 )} 
a-2i(p, k) 

(fc2 l p2) i{{k-pf + k 2 {k-p) 2 -ik 2 p 2 k-p 
+2k 4 k-p + k 4 p 2 -2k 6 }m 2 /k 2 

+(k ■ pf + (2k 2 - p 2 )(k ■ pf + k 2 p 2 k ■ p - 2k 4 k ■ p - k 2 p 4 
M£- 1) {(k ■ pf +P 2 (k ■ pf-3k 2 p 2 k ■ p 

+2k 4 k-p+k 2 p 4 -2k 4 p 2 } 

«23(P, k) 

q 2 (m 2 + k ■ p) 

+ (£ - l)(g 2 m 2 + (k-pf - (k 2 +p 2 )(k ■ p) + k 2 p 2 ) 



X 



q 2 k ■ pm 4 + 2(k 2 + p 2 )A 2 m 2 



-k 2 p 2 {2(k ■ pf + (k 2 + p 2 )k ■ p - Ak 2 p 2 } 



' 2 



■ m 



{-2(k ■ pf + k 4 + p 4 } 
+2 {(k 2 + p 2 )(k ■ pf + (k 2 - p 2 f k-p- k 2 p 2 (k 2 + p 2 )} to 2 
+ \ {-4(fc • pf + (k 2 + p 2 f(k ■ pf + k 2 p 2 (k 2 - p 2 ) 2 } 



+ 



(£-1) 



q 2 {-2{k ■ pf + k 4 + p 4 } m 8 



+2 



{(fc 2 + p 2 )[-2(k ■ pf + (k 2 +p 2 )(k- pf + (k 4 + p 4 )k ■ p] 
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a-32(k,p) 



a i2 {k,p) 

a 4 3(k,p) 

au{k,p) 
a 45 (k,p) 

a 4Q (k,p) 



-fcV(3fc 4 + 3/-2fcV)}™ 6 

-^q 2 A 2 (k 2 -p 2 ) 2 m 4 - 2{(k 2 + p 2 )(k 4 + p 4 - 4fcV)(fc • pf 
-k 2 p 2 (k 4 + p 4 - 6fcV)(fc ' pf 
+2k 4 p 4 (k 2 + P 2 )k ■ p - k 4 p 4 (k 2 + p 2 ) 2 }m 2 

-\k 2 p 2 q 2 {(k 4 +p 4 - Qk 2 p 2 ){k-p) 2 + k 2 p 2 (k 2 +p 2 ) 2 ) 

<X3l(f>, k) 

i{{k- pf~k 2 {k ■ p) 2 -Zk 2 p 2 k-p+2k 4 k-p-k 4 p 2 + 2k 6 } m 2 /k 2 
+ {(i-2){{k-pf-{2k 2 -p 2 ){k-p) 2 







a 34 (k,p) = 


= CL33(k,p) 


d3h{k,p) = 


- -(fc 4 +p 


a3e(k,p) = 


= -mA 2 |/ 




+ «-D 



„2 i „2\ , ol,2^2 



'221 4 

1 m 



a 4 i(fc,p) = -771 (C-l) 



< Z 2 {(fc 2 +p 2 )fc-p+2fcV} 
+2(fc 2 + p 2 ) 2 A 2 m 2 - k 2 p 2 (k + p) 2 { (k 2 + p 2 )k ■ p - 2k 2 p 2 } 
(k 2 -p 2 ) 



2x 



q 4 m 6 + 3g 2 {-(k 2 +p 2 )k ■ p+2k 2 p 2 } m 4 

+ {(*;■ p) 2 [4(fc • p) 2 - 3fc 4 - 3p 4 - 26fcV] 
+fcV[24(fc 2 + p 2 )fc • p - 3fc 4 - 3p 4 - 14* Y)]}^ 
+ ^{(fc 2 +p 2 )(fc-rt 3 + 2fcV(fc-p) 2 

-3fcV(fc 2 +p 2 )fc-p + 2fcV} 

-04l(p,fc) 

[(fc 2 + fc-p)(fc-p) 2 + fc 2 (2fc 2 -3p 2 )fc-p+fc 4 (p 2 -2fc 2 )] 
-a 43 (p, fc) 



= (e-l)(fc 2 -pV 



A 2 



mg - 1) (fc 2 - p 2 ) — \q 2 k ■ pm 2 + 2(k 2 + p 2 ) (k ■ p) 2 
X 
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-2k 2 p z k-p-k 2 p 2 (k 2 +p 2 " 
1 4x 



a 51 (k,p) 



a^{k,p) 

a 54 (k,p) 
a 55 (k,p) 

a 56 (k,p) 



a&2{k,p) 
a&z{k,p) 



a M (k,p) 
a 65 (k,p) 



+6q 2 {2k 2 p 2 - (k 2 + p 2 )k ■ p) m 4 - 6k 2 p 2 q 4 m 2 

-q 2 {(k 2 -p 2 f(k ■ p) 2 + 2k 2 p 2 {k 2 +p 2 )k ■ p-k 2 p 2 (k 2 +p 2 ) 2 } 

«5i (p, k) 

(L^l[(k-p) 2 + 2k 2 k-p-k 2 (2k 2 +p 2 )} 

«53(P, k) 

(e-i)g 2 

-m(£-l) — [q 2 (k 2 +p 2 )m 2 + 2(k 4 +p 4 )k-p-2k 2 p 2 (k 2 +p 2 )] 



a e i(k,p) = -T]i 



(k 2 -p 2 ) 



q 2 m 4 - 2 { (k ■ p) 2 - (k 2 + p 2 )k ■ p + k 2 p 2 } to 2 



+\{{k-vf + k 2 p 2 } 



+ 



g 2 (C-l) 

X 



q 2 m 8 + 2 { (k ■ p) 2 + (k 2 + p 2 )k ■ p - 3k V } 



in 



--q 2 A 2 m 4 -2{(k 2 + p 2 )(k ■ pf - k 2 p 2 (k ■ p) 2 - k 4 p 4 } m 2 

-]-k 2 p 2 q 2 {{k-p) 2 + k 2 p 2 } 
-a 61 (p, k) 

- £ {{k 2 + k ■ p)(k ■ p) 2 + k 2 (2k 2 - 3p 2 )k ■ p 
~k 4 (2k 2 -p 2 )}m 2 /k 2 

+ {^-2){{2k 2 -p 2 + k-p){k-p) 2 -k 2 {2k 2 -p 2 )k-p~k 2 p 4 } 
-a&3(p, k) 

-q 2 (k 2 -p 2 ) [Cm 2 -(C-2)k-p] 



a 66 (k,p) = -m(k 2 -p 2 )A 2 
— k 2 p 2 q 2 k ■ p) 



k-p+ ^LJl (q 2 k ■ pm 4 + 2{k 2 + p 2 )A 2 m 2 
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a n (k,p) 



qi(g-l) 
4x 



- 2q 6 m 6 - 6g 4 {(k 2 + p 2 )(k ■ p) - k 2 p 2 } 



2 2\ 4 

1 m 



a72(k,p) 

a73(k,p) 
a 74 (k,p) 
a 75 (/c,p) 



-3q 2 {{{k-p) 2 + k 2 p 2 ){k A +p i + Qk 2 p 2 )~Sk 2 p 2 {k 2 +p 2 )k-p} : 
+9 2 {(fc 2 - P 2 ) 2 (fc • P? + ik 2 p 2 (k 2 + p 2 )(k ■ pf 

-k 2 p 2 {ik A + 3p 4 + Wk 2 p 2 )k ■ p + Ak A p i {k 2 + p 2 )} 

on(p, fc) 

(£ - l) (fc2 ~ 2 fc " P) ((fc ■ p) 2 + 4fc 2 fc • p - 2fc 4 - 3fcV) 



fc 2 



«73(P, fe) 



a 7& {k,p) = m(£-l)- 



<? 2 <^ (k 2 +p 2 )k-p-2k z p z }m 



2„2 I 2 



+2(fc 4 +p 4 )(fc • p) 2 -4fcV(fc +P )fc ' p-fcV(fc +P -6/cV) 



ffl8l(fe,p) 

as2(k,p) 
a%z(k,p) 
a S i(k,p) 
as5(k,p) 
a S6 (k,p) 



= -vi- 



-q 2 (m 2 + k ■ p) 



081 (p, fc) 

2(£ + 2)fc- g 

OS3(P, fc) 

-(C + 2)g 2 
. 



(5.163) 



An important point to note is that these coefficients do not contain any 
trigonometric function, as it has been extracted out for raising the Tj to a non- 
perturbative status. The Tj have the required symmetry under the exchange of 
vectors k and p. All the Tj are symmetric except T4 and T6 which are antisym- 
metric. Note that the form in which we write the transverse vertex makes it 
clear that each term in all the Tj is either proportional to al(l) or a/(k 2 p 2 ). We 
shall see that this form provides us with a natural scheme to arrive at its simple 
non-perturbative extension. 

A few comments in comparison with the work by Davydychev et. al. |67| . 
are as follows: (i) None of the Ti we have calculated has kinematic singularity 
when k 2 — > p 2 . This clearly suggests that the choice of the basis {T^} suggested 
by Kizilersu et. al. is preferred over the one of Ball and Chiu (in QED3 as well) 
used by Davydychev et. al. 67 . In particular our T4 and T7 are independent 
of kinematic singularities, (ii) In three dimensions, their factorization of the 
common constant factor in Eq. (E.l) is singular. However, as the divergences 
completely cancel out, we find our expressions more suitable for writing the 
transverse vertex in three dimensions, (iii) With the way we express Jo, all the 
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Tj are written in terms of basic functions of k and p and a single trigonometric 
function of the form 1(1). This form plays a key role in enabling us to make 
an easy transition to the possible non-perturbative str ucture of the vertex, as 
explained in the next section. Moreover, with the given form of Jo, a direct 
comparison can be made with the massless case. 



5.3 Non-perturbative Form of the Vertex (Ef- 
fective Transverse WTI) 

5.3.1 On the Gauge Parameter Dependence of the Vertex 

Let us first look at the Tj in the simplified massless case, with the notation 
k = \/ — k 2 , p = \] — p 2 and q = \/—q 2 |571 , 



T-2 



T 3 = 



an 



4 kp(k + p){k + p + q) 2 



i + K-i) 



2k + 2p + q 



(5.164) 



[4kp + 3kq + 3pq + 2q 2 



8 kpq(k + p + q) 2 
+ - 1) (2k 2 + 2p 2 + kq + pq)] , 



1~6 = 



k — p 



a7r(2 - Q 

8 kp(k+p + q) 2 ' 

a7r(2 + £) 1 



2 kp(k+p + q) ' 
It is interesting to note that the existence of the factor 

k-p (\ V 



(5.165) 
(5.166) 

(5.167) 



kp 



k 



P 



in eq. H5.166|l puts tq on a different footing as compared to the rest of the Ti 
The reason is that in the massless limit, the Fermion Propagator is simply 



F(p) 



= 1 



7ra£ 1 
~ p ' 



implying 



1 



F(k) F(p) 



Therefore, the relation of tq with the Fermion Propagator of the type [l/F(k 2 ) 
1/F(p 2 )} seems to arise rather naturally : 



T6 



2-e 



2£ {k+p + q) 2 



F(k) F(p) 



(5.168) 
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as noticed first by Curtis and Pennington [37j (note however that their coefficient 
is not the same). In the rest of the Tj, the factor 1/k — 1/p docs not arise. 
However, one could introduce it by hand to arrive at the following expressions : 



T-2 



'£ (k 2 - p 2 )(k + p + q) 2 I ! " i; " 11 



2k + 2p + q 

q 



T3 



T» = 



F(k) F(p) 

1 [Akp + 3kq + 3pq + 2q 2 + (£ - 1) (2k 2 + 2p 2 + kq + pq)] 



(5.169) 



2£ 



1 



F(k) F(p) 
2(2 + 



1 



q(k - p)(k +p + q) 2 



1 



£ (k-p)(k+p + q) 



1 



F(k) F(p) 



(5.170) 
(5.171) 



Equations. I|5. 16815. 171|l represent a non-perturbative vertex which is in agree- 
ment with its complete one-loop expansion. This vertex has been constructed 
in accordance with the form advocated, e.g., in [7|E|021- There are a couple 
of important points which need to be discussed here : 

• There is an explicit dependence on the gauge parameter, £. A widespread 
belief has been that the gauge dependence of the vertex should solely 
arise through functions F(k 2 ) and F(p 2 ), and there should be no explicit 
appearance of the gauge parameter £. Such a belief has been expressed 
(or is reflected) in various works to date, e.g., [7[ El EZ1 EH EH1 ■ Here we 
show that at least in massless QED3, such a construction is not possible. 

• One of the main reasons that the transverse vertex was believed to be 
proportional to the factor 



1 



1 



F(k) F(p) 

was the assumption that the transverse vertex vanishes in the Landau 
gauge, based upon perturbative results. A complete one-loop calcula- 
tion reveals that the transverse vertex does not vanish in the Landau 
gauge. Moreover, an explicit presence of the gauge parameter in the non- 
perturbative form of the vertex tells us that the presence of the factor 
[1/F(k) — 1/F(p)\ is no longer a guarantee that the transverse vertex 
vanishes in the Landau gauge. 

We now show that the explicit dependence of the vertex on the gauge pa- 
rameter £ is unavoidable in massless QED3. We notice that at the one loop 
level, each of the r% can be written in the following form : 



n(k,p,q) = aai(k,p,q) + a£bi(k,p,q) 
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On the other hand, Eq. I|5.5[) yields the following form for F 

1 



F(p) 



1 + Ci(p) 



If we want to write the non-perturbative form of the Tj in terms of 1/F(p) and 
l/F(k) alone and we do not expect explicit presence of a, the only way to get 
rid of £ dependence is to have 

b 2 T% + b 3 T 3 " + b 6 Tg + b & T g " = . 

It is not possible as Tjf form a linearly independent set of basis vectors. There- 
fore, any construction of the 3-point vertex will surely have an explicit depen- 
dence on the gauge parameter. Owing to these reasons, we realize that to 
demand the transverse vertex to be proportional to [1/F(k) — 1/F(p)] is artifi- 
cial (apart from tq) and is not required. Therefore, we do not pursue this line 
of action anymore. In the next section, we move on to construct the vertex for 
the massive case inspired from our perturbative results. 



5.3.2 Non-perturbative Vertex (Effective Transverse WTI) 



As pointed out in the previous section, each term in all the Ti is either propor- 
tional to the trigonometric function al(l) or a/(k 2 p 2 ). On the other hand, the 
perturbative expressions for A4(p) and F(p), Eqs. 1(5.5(1 . permit us to write : 



F(k) 



where 



and 



2(2+£)RT(Z) 



F(p) 
A(p) 

M(l) 
F(l) 



fcfaf [k 2 A(p)- P *A(k)}, 



{to — (to 2 + p 2 )I(p)} 



(5.172) 



(5.173) 



F(l) 



£(m 2 + l 2 ) 
21 2 



alii) . (5.174) 



In the massless limit, Eq. (|5.172|) simply reduces to 



F(k) F(p) 



in the Euclidean space, as expected. It was in fact an analogous massless ex- 
pression in the limit when k > > p that inspired Curtis and Pennington, |37| , to 
propose their famous vertex in QED4. Here, we are extending the reasoning to 
all the momentum regimes in the massive QED3. Fortunate simultaneouss oc- 
currence of the factor a/ (k 2 p 2 ) in all the 8 Eqs. (|5.162() and Eq. I|5. 172(1 . and the 
presence of the same trigonometric factor 1(1) in the expressions for the vertex as 
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well as the propagator, one naturally arrives at the following non-perturbative 
form of n : 



2a i6 (fc,p) 



( MM 

2(£+2)El(l 3 -) \ F(lj) 



-m - 1- 



1 



F(lj) 



+ 



£[k 2 A{p)-p 2 A(k)] 



F(k) F(p) 



(5.175) 



By construction, in the weak coupling regime, this non-perturbative form of 
the transverse vertex reduces to its corresponding Feynman expansion at the 
one loop level in an arbitrary covariant gauge and in all momentum regimes. 
Note that we have managed to write the transverse vertex solely as a func- 
tion of the fermion propagator. Therefore, effectively, we have a WTI for this 
part of the vertex. We would like to emphasize that this is not a unique non- 
perturbative construction. However, it is probably the most natural and the 
simplest. A two loop calculation similar to the one presented in this chapter, 
and the LKF transformation for the vertex should serve as tests of Eq. I|5. 175(1 or 
guides for improvement towards the hunt for the exact non-perturbative vertex. 
On practical side, the use of our perturbation theory motivated vertex in studies 
addressing important issues such as Dynamical Mass Generation for fundamen- 
tal fermions should lead to more reliable results, attempting to preserve key fe 
atures of gauge field theories, e.g., gauge independence of physical observables. 
A computational difficulty to use the above vertex in such calculations could 
arise as the unknown functions F and M. depend on the angle between k and 
p. This would make it impossible to carry out angular integration analytically 
in the SDE for the Fermion Propagator. This problem can be circumvented 
by defining an effective vertex which shifts the angular dependence from the 
unknown functions F and M to the known basic functions of k and p. This can 
be done by re- writing the perturbative results, Eq. (|5. 162(1 . as follows : 



where 



n(k,p) 

bn(k,p) 
b i2 (k,p) 



b a (k,p)I(k) +bi 2 (k,p)I(p) 



ai6(k,p) 
k 2 p 2 

an{k,p) -—- + a i3 (k,p) + - a l5 (k,p)—— 
1(h) ^ 1(h) 

ia 2 ) . . i „ A(iT) 

a i2 (k,p) ——t + a i4 ,{k,p) + - a l5 {k,p)—— 

1(14,) * 1(H) 



(5.176) 

(5.177) 
(5.178) 



This form can now be raised to a non-perturbative level exactly as before, with 
the only difference that the functions F and M. are independent of the angle 
between the momenta k and p : 



f*( 2b ij (k,p) K ? j 
9l \U\^m 2 + n 2 ) 



M(Kj) 
2(£ + 2)/ £ ?J(« i ) \F(kj) 
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- 1- 



1 



where k± = k and K2 = P 



2a i6 (k,p) 



£[k 2 A(p) - P 2 A(k)} 



F(k) F(p) 



(5.179) 



5.4 Comments on the d— dimensional Case 

Following Davydychev et. al, (HZ), two of the Master Integrals needed to cal- 
culate are : 



(5.180) 
(5.181) 



1(^1,^2,^3) = J 



[(p - w) 2 — m 2 ]" 1 [(k — w) 2 — m 2 ]^ [w 2 ]^ 
d d w 

[(p - w) 2 ] Ul [(k - w) 2 }^ [w 2 - m 2 }^ 



In terms of these integrals, one can write out the Fermion Propagator and 
fermion-boson vertex in the context of QED in arbitrary dimensions. 

5.4.1 Fermion Propagator 

Fermion Propagator in the d-dimensional case at the one-loop order can be 
written as : 



1 



e 2 (d-2)£ r , 2 , 2 



F(p) i{2n) d 2p 



, [(p 2 +m 2 )I(0,l,l)-Z(0,0,l)] , 



M(p) e 2 
F(p) i{2n)' 



Md- 1-01(0,1,1) 



(5.182) 



The integral 1(0, 1, 1) is the <i-dimensional counterpart of our I(y) function in 
the corresponding equation in three dimensions. 

5.4.2 The Transverse Vertex 

The corresponding d-dimensional r's of the transverse vertex can be expressed 
as, IHZI : 



i{2ir) a I 

+ g 2 /(l, 1, 0)] + t i>2 (l(0, 1, 1) + 1(1, 0, 1) - 2/(1, 1, 0)) 

+ u, 3 (i(o, 1, 1) + 1(1, 0, 1) - 2 X(0 ' °; 1} 



+ ^4(1(0, 1,1) +1(1, 0,l))+t ( , 8 



T(0,1,1)-J(1,0,1) 

fc 2 — p 2 



(5.183) 
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where the ti t j, j = 1 . . .5 are basic functions of k,p and q, listed in |S7]. A 
straightforward procedure to obtain a non-perturbative vertex in this case is a 
natural thing to ask for. However, a dificulty arises since the function 7(1, 1, 1) 
in general cannot be rewritten as a linear combination of the integrals which 
appear in the Fermion Propagator in an obvious way. Note that in threee 
dimensions 7(1,1,1) has the same structure as the integrals which define the 
Fermion Propagator, though with more complicated arguments. Therefore, we 
have not yet been able to write a simple non-perturbative extension of the 
transverse vertex in arbitrary dimensions. 

5.5 Comments on the Two Loop Case 

As we have mentioned earlier in this chapter, the three point vertex must reduce 
to its all order Feynman expansion in the weak coupling regime. What we have 
achieved in our work is the realization of this statement at the one loop level. To 
be able to do so beyond this order, we must know the perturbative vertex at the 
higher loops to have a guide for the construction of its non-perturbative coun- 
terpart. A good amount of progress has already been made in the calculation 
of necessary integrals involved. Most of these attempts surge from the precision 
analysis in particle physics phenomenology at the next-to-next-to-leading order 
calculations in perturbation theory. However, the scattering (see for example 
or decay processes (see for example |HH] ) do not in general require the 
calculation of vertex- type integrals with off-shell fermions and bosons. A com- 
plete two-loop calculation of the transverse vertex in arbitrary covariant gauge 
in arbitrary dimensions, to the best of out knowledge, still d oes not exist for 
off-shell particles in the legs. 

Using Dimensional regularization with d = 4 — e, the integrals appearing in 
the calculation of two-loop corrections take the generic form |S7] : 



where Pi are the propagators, depending on W\, W2 and the external momenta 
Pi . . . p n while Si are scalar products of a loop momentum with an external 
momentum or of the two loop momenta. The integral is specified by the powers 
rm, (mi > 1) of all propagators (Pi, . . . , P t ) and by the section (Si, . . . , S q ) of 
scalar products and their powers (ni . . . n q ), (n, > 0). 

Progress in the calculation of two-loop corrections to four-point amplitudes 
was based on an efficient procedure to reduce the large number of different 
scalar integrals to a very limited number of so-called Master Integrals. The 
scalar integrals are related among each other by various identities. One class of 
identities are the so-called integration by parts (IBP) identities : 




(5.184) 




(5.185) 
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where V is any combination of propagators, scalar products and loop momen- 
tum vectors. V can be a vector or tensor of any rank. IBP identities are a 
generalization of (|5.92l) and they follow from the fact that the integral over the 
total derivative with respect to any loop momentum vanishes. 

Another class of identities is obtained from the fact that all integrals under 
consideration are Lorentz scalars, which are invariant under Lorentz transfor- 
mations of the external momenta |S8) ■ These Lorentz invariance (LI) identities 
are obtained from : 

Master Integrals can be analytically expressed in series of the Dimensional 
Regulator parameter e using Feynman or Schwinger parametrizations. However, 
sometimes using these parametric representations does not leave an integral that 
can be solved easily. Different techniques arose due to this situation, among 
others : 

• The Mcllin-Barnes method is based on the representation for a sum to 
some power, as a contour integral over a complex variable and the integra- 
tion is then performed on straight contour lines parallel to the imaginary 
axis. After closing the contour, the result is the sum of all enclosed residues 
that might be expressed as a hypergeometic series. For the two-loop boxes 
we have [55] , 

• The Negative Dimensions technique consists in rewriting the integral over 
the parameters by introducing new ones through a multinomial expansion. 
Many conditions have to be satisfied among the parameters, which leads 
to the restriction that d must be a negative integer. Some results arc 90 . 

• Double-Integral representation, that consist on splitting the loop space 
into the parallel and perpendicular subspaces to the external momentum. 
Then, using cylindrical instead of spherical coordinates, one can calculate 
the real and the imaginary parts of the integrals separately, performing 
contour integrals. For some results, see |91| . 

• Strategy of Expansion by regions , in which, instead of the integration 
over the whole space of loop momenta, the integration is performed only 
over some specific regions. The strategy consists in considering various 
regions of the loop momenta and expand, in every region, the integrand 
in Taylor series with respect to the parameters that are considered small 
in the given region; integrate the integrand expanded in every region and; 
put to zero any scaleless integral. 

• Dispersion relations are used to calculate self-energies for the small exter- 
nal momentum behavior |93j . 

• Numerical strategies are used, where any loop integral is stripped analy- 
tically of its IR singularities so that the finite integrals can be performed 
numerically. |94| . 
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• The analytic evaluation of MI can also be carried out without explicit 
integration over loop momenta by deriving differential equations for MI in 
internal propagator masses or external momenta. The equations can be 
solved with appropriate boundary conditions |95j . 

5.5.1 Two-Loops Fermion Propagator 

The two-loops contributions to the fermion propagator are depicted in the fol- 
lowing Diagram : 



Diagram 9 : Two-Loop Corrections to the Fermion Propagator. 

Using the Double-Integral representation method, the Fermion Propagator has 
been calculated in [5T] for arbitrary masses. Scalar self-energies have been cal- 
culated using dispersion relations ||3| ■ Using the Pinch Technique [HHI > off-shell 
self-energies have been calculated, in terms of Harmonic Polylogarithms. Also, 
an asymptotic expansion has been used to calculate the on-shell master dia- 
gram for the two-loop propagator This expansion is taken in the limit of 
the external momenta, and therefore the fermion masses, tend to infinity. On- 
shell self-energies master integrals with one mass have been calculated in |98| . 
In the light-cone gauge, the fermion propagator has been calculated in |9T)] . 
Following the algorithm developed in |100| . these corrections have been calcu- 
lated in an arbitrary covariant gauge for off-shell massive fermions in arbitrary 
dimensions jlOlj . 
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5.5.2 Two-Loops Vertex 

The two-loops corrections for the vertex are shown in Diagram (10) 
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Diagram 10 : Two-Loop Corrections to the Vertex. 

Some partial calculations to the full vertex using the double integral repre- 
sentation, planar vertex in four dimensions has been presented in [91| . In Ref. 
|l(J2j a recent on-shell calculation for the vertex has been done at arbitrary 
momentum transfer in arbitrary dimensions. The planar vertex ((a) in Figure 
5) with essential on shell singularities has been calculated in arbitrary dimen- 
sions |1U3| applying the strategy of expansion by regions. The calculation of a 
the vertex function at zero momentum transfere and with space-like external 
momentum |1(J4| . Vertex with two legs on-shell in arbitrary dimensions has 
been calculated in |l(J5j using the Negative Dimensions technique. Numerically 
speaking, the massless vertex has been calculated in the Landau gauge |106| . 
For non-planar vertex diagrams ((b) in the corresponding figure), a numerical 
method is propossed in jl()7| . In the Sudakov limit, the Vertex has been cal- 
culated using the Mellin-Barnes method in arbitrary dimensions [108) . A full 
off shell calculation of the two- loops vertex has not been carried out yet. A 
programme towards the construction of the full two-loop vertex has to take into 
account the transverse part of it, so that one can check of modify the corres- 
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ponding basis {T/ 4 }. With this basis one can look for a nonpcrturbative vertex 
in arbitrary dimensions, or, at least, in our three dimensional case. 

Once we are able to construct a three-point vertex which agrees with its 
perturbative expansion in the weak coupling regime to all orders, we shall rest 
assured that all the gauge identies are satisfied at every order of perturbation 
theory. We can then hope that the physical observables associated with non- 
perturbative physics shall automatically be gauge independent. 



Chapter 6 

Discussion and Conclusions 



In this thesis we have focused our atention in incorporating gauge identities 
(Ward-Green- Takahashi identities and Landau-Khalatnikov-Fradkin transfor- 
mations) to the non perturbative study of the Schwingcr-Dyson Equations in 
QED3. We know that Perturbation Theory is the only scheme where in each 
order of approximation, such identities are satisfied, and also it verified that, 
indeed, physical observable quantities are independent of the gauge parameter. 
However, in a non perturbative analysis of the Schwinger-Dyson Equations it 
has been impossible the conjunction of these two facts. 

A trace back for the origin of the gauge dependence points out to the fol- 
lowing problems : 

1. Violaton of the Ward-Green- Takahashi Identity. 

2. Incorrect gauge behavior of the Green's functions. 

3. Use of wrong technical tools, like gauge dependent regularization schemes. 

In order to avoid that the technical details overlap the pure theoretical issues 
of the study, we choose QED3, because in a theory that lacks of ultraviolet 
divergences. Therefore, the source for the gauge dependence should be only due 
to the first two points above mentioned. 

In Chapter 3 we considered the simplest scenario, and we calculated the 
mass and the chiral condensate by solving the Schwinger-Dyson Equation for 
the Fermion Propagator. We explicitly observed that these quantities depend 
upon the gauge parameter. Since an incorrect regularization cannot induce such 
dependence, the only room for improvement is in the three-point vertex ansatz. 

An ansatz should be such that : 

• It satisfies the Ward-Green- Takahashi Identity, 

• The Fermion Propagator obtained from this ansatz and the vertex itself 
should satisfy their corresponding Landau-Khalatnikov-Fradkin transfor- 
mation. 
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The incorporation of the Ward-Green- Takahashi Identity is straightforward. 
However, as we did in Chapter 4, even in the simplest case, the Landau- 
Khalatnikov-Fradkin transformation for the Fermion Propagator leads to a 
gauge behavior extremely complex for this Green's function, and in spite of 
this, we were able to extact valuable information about its structure beyond the 
tree level. The search for a non perturbative vertex which leads to this gauge 
bahavior for the Fermion Propagator is not at all a trivial task. It is here that 
Perturbation Theory plays its role. If we can put forward to a non perturba- 
tive ansatz which reduces to its Feynman expansion to all orders in the weak 
coupling regime, we can hope that it automatically incorporates the Landau- 
Khalatnikov-Fradkin transformations not only for the Fermion Propagator, but 
for the Vertex itself. 

Therefore, in Chapter 5 we calculated the Fcrmion-Boson Vertex at one-loop 
level outside of the quiral phase of QED3, previously studied. We decompose 
this vertex into its longitudinal and transverse parts. The former guarantees 
the validity of the Ward-Green- Takahashi Identity, since it is related to the 
Fermion Propagator calculated at the same order of approximation, while the 
later ensures that the Landau-Khalatnikov-Fradkin transformations are indeed 
satisfied. Taking advantage of the lack of additional constraints, we exploited 
the perturbative form of the Fermion Propagator in such a way that we could 
write an effective Ward-Green- Takahashi Identity for the transverse vertex. This 
identity allows us to have the first insight on the non perturbative structure of 
the Fermion-Boson interaction, since the transverse vertex, related in this way 
to the Fermion Propagator does not depend explicitly on the electromagnetic 
coupling. We write, afterwards, the transverse vertex in a convenient form that 
will allow, to corroborate its non perturbative structure in a future work at the 
two-loop level. 

The systematic construction of the three point vertex guarantees that : 

• The Ward-Green- Takahash identity, which relates the Fermion Propagator 
to the Fcrmion-Boson vertex, is satisfied non-perturbatively. 

• Both the vertex and the resulting Fermion Propagator satisfy Landau- 
Khaltnikov-Fradkin transformations to 0(ct) and 0(a 2 ) respectively. 

• The vertex does not contain any kinematic singularities when k 2 — > p 2 . 

• Most importantly, the vertex reduces to its correct perturbative expansion 
at 0(a) in the weak coupling regime in an arbitrary covariant gauge and 
for all momentum regimes. 

• The vertex has correct symmetry under the parity charge congugation 
and time reversal operations. 

This will allow that the technique as well as the reasoning used in our cons- 
truction could be applied in more complicated theories like QCD in a suitable 
and realistic way. Finally, it is also possible to incorporate our ideas in some 
alternative scheme to the Standard Model which pretends to give a solution to 
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the problem of the origin of the masses, 
whishes to do so. 
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